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SUMMARY 


Procedures are reco mme nded for solving tie equations of equilibrium 
of reinforced panels and isolated fuselage rings as represented ly tie 
external loads and tie operations talle estallisled according to 
Southwell’s metlod- From tie solution of tlese equations tie stress 
distribution can le easily determined. Tie recommendations are lased on 
tie experience of tie past 1 years In applying numerical procedures to 
monocoque stress analysis at tie Polyteclnic Institute of Brooklyn 
Aeronautical laboratories . Tie metlod of systematic relaxations, tie 
matrix calculus metlod, and several otler metlods applicable in special 
cases are discussed. 

Definite recommendations are made for obtaining tie solution of 
reinforced -panel problems -which are generally designated as shear lag 
problems. The procedures recommended are demonstrated in the analysis 
of a number of panels, several of which were discussed in previous 
PIBAL reports, whereas others are shown for the first time. 

In tie case of fuselage rings it Is not possible to make definite 
recommendations for tie solution of the equilibrium equations for all 
rings and loadings • However, suggestions based cm. the latest experience 
are made and demonstrated on several rings. 


INTRODUCTION 


The application of the indirect metlods of Hardy Cross (reference l) 
and R. Y . Southwell (reference 2) to tie analysis of monocoque structures 
las been shown In a series of investigations (references 3 to 8) carried 
out at tie Polyteclnic Institute of Brooklyn Aeronautical Laboratories. 
Tlese indirect methods are likely to lead to solutions of problems in 
stress analysis that are intractable by direct analytical metlods because 
the structure is tapered, las large cut-outs, its reinforcing elements 
are distributed irregularly, or tie like. 

The distorted shape corresponding to equilibrium under the applied 
loads is determined first in tie indirect metlods. From It the stresses, 
forces, and moments required can be calculated without difficulty. This 
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approach is Justified "by the comparative ease with which the stresses in 
a coup lei structure can he determined for an individual displacement of- - 
one point and with which the final distorted shape of a coup lex structure 
can he represented by a summation of such individual displacements. 

The complete structure is considered to he composed of appropriate 
elements and its degrees of freedom are the displacements of the several 
reference points on the boundary of each elements Each of these points 
is displaced in turn and the reactions at the reference points caused hy 
the displacement are listed. If hy suitable displacements - of nl 1 points 
the reaction forces and moments are made equal and opposite to the 
external loads at each point, the Whole structure is in equilibrium and 
its distorted shape is determined. 

In applying' the indirect methods to monocoque structures the termi- 
nology of Southwell (reference 2 ) has been retained. Thus, the elements 
which compose the complete structure are "'units " and the determination 
of the forces and moments due to" a displacement of a boundary point of 
such units is termed the 'Wit problem. " The magnitudes of these forces 
and moments are given by "influence coefficients. " The complete effect 
of a displacement is given in an "operations table, ” and the step-by- 
step process, which can be employed to determine the equilibria 
distorted shape, is called the ’iaethod of systematic relaxations. ” At 
each step of this process forces and moments referred to as "residuals " 
remain unbalanced at - each point in the structure. A running account of - 
the residuals and of the displacements or "operations " undertaken is 
kept in the '^relaxation table. " 

The operations table along with the external forces constitutes a 
system of linear equations, which are equal in number to the degrees of 
freedom of the structure and which have as variables the displacements. 
Each equation represents the condition of - equilibrium for the force or 
moment associated with one degree of freedom. When the method of system- 
atic relaxations is applied an approximate solution to this system of 
equations and accordingly an approximate - equilibrium state - of the struc- 
ture are found. 

The indirect method of analysis Just— outlined has been applied at 
PIBAL to the reinforced -panel and ring components of a monocoque struc- 
ture as well as to complete - circular cylinders with and without cut-outs. 
In references 3 and 4 the stress distribution in the sheet - and stringers 
of a reinforced panel was determined under loads applied parallel to the 
stringers. Fuselage rings with and without internal bracing elements 
were investigated in reference 5. The determination of- the influence 
coefficients for the ring unit problem was found to involve considerable 
computational work and therefore appropriate graphs and tables are 
given in reference 6 to facilitate their calculation. In references 7 
and 8 the elements, namely, the reinforced panel and the ring, are 
combined into a circular cylinder and the stress distribution in the 
cylinder was investigated for the case when the loading is a pure bending 
moment • 
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In the application of the indirect -stress -analysis methods to the 
problems mentioned the major obstacle has been to find an appro ximat e 
solution of the system of equations -with a reasonable expenditure of 
effort • In each problem it has been readily possible to establish 
satisfactory units and to combine them to represent the con® lex struc- 
ture. During the past h years considerable experience has been gained 
at PIBAL in overcoming this obstacle to the wider application of numerical- 
procedures in the analysis of monocoque structures. On the basis of this 
experience same recommendations can be made as to the most expeditious 
method of solving reinforced -panel and fuselage -ring problems after the 
operations table has been established as described in references 3 to 5* 

In many problems solution of the set of linear equations by means of 
matrix algebra was found easier and less time consu m i n g than the 
determination of the displacements by systematic relaxations • In other 
canes special methods, such as the growing -unit method, proved to be 
most expeditious • 

It Is assumed that the reader Is familiar with the terminology of 
Southwell's relaxation method and with the solution of the unit problem 
as well as the establishment of the operations table for both the 
reinforced -panel and fuselage -ring problems. Couplets details of these 
are given in references 3 to 6. 

ThlB work, carried out at the Polytechnic Institute of Brooklyn, was 
sponsored by and conducted with the financial assistance of the National 
Advisory Committee for Aeronautics. Arnold 0. Ostrand contributed the 
growing-unit method for reinforced panels. The authors also wish to 
acknowledge their indebtedness to the following members of the staff of 
the Polytechnic Institute of Brooklyn: Professors George B. Hoadley and 

William Maclean of the Department of Electrical Engineering for their 
work on the electric analogue, Burton Erickson for carrying out the 
major portion of the computations, and Bruno A. Boley for his editorial 
advice • 


SYMBOLS 


A cross-sectional area of stringer and effective sheet 

A - Q points . on a ring or a reinforced panelj group operations 

A* effective shear area of ring section 

a distance between adjacent longitudinal stringers 

b distance between adjacent transverse stringers 

C 


electrical conductance 
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Young r s modulus of elasticity 

tensile force in stringer} applied external load 
shear modulus of elasticity 
horizontal direction 

moment of inertia of cross section} electrical current 
group operations 

length of straight-bar or length of arc of curved "bar 

bending moment 

moment acting on a Joint 

shear flow 

radial force acting at a joint} electrical resistance 
tangential force acting oma Joint 
sheet thickness 

displacement of a Joint in tangential direction 
electrical potential} vertical direction 

displacement of a Joint in radial direction} displacement of' 
a joint in vertical direction 

vertical "block: displacement 
rotation of a joint 

magnitude of group operation to be determined 

rectangular coordinates 

force in y-axis direction 

angle subtended by ring segment 

sect ion -length parameter (AL^/l) 

ratio of effective shear area to tension area (A*/A) 
summation 
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REINFORCED PANELS 
Introduction 


In this section plane and slightly curved reinforced panels are 
discussed -when the loads are applied in the plane of the flat panels or 
tangentially to the surface of the slightly curved panels. 

In most airplane structures there is a predominant direction in 
■which the major forces act and In -which the major reinforcing elements 
lie. When the panel is symmetric and symmetrically loaded experience 
has shown that it suffices to consider displacements and force equi- 
librium in the predominant direction only. Even when the structure or 
the loads are nonsymmetric, the displacements and forces in the trans- 
verse direction are usually of secondary importance hut they may he 
considered in a more refined analysis. 

In references 3 and- 4 numerical procedures for the determination of 
the stress distribution in reinforced panels subjected to axial stringer 
loads are developed and demonstrated on several flat and curved panels 
with and without cut-outs. The results obtained by means of these 
procedures are in good agreement with those of tests. 

Solution of the system of equations represented by the operations 
table and the external forces can be found by several methods, five of 
which are described herein. The various conditions of loading and 
structure which suggest the use of one method rather than another are 
discussed. 


Relaxation Method 

For most reinforced-panel problems the relaxation method of 
solution is the most suitable. Simple group and blocfc. operations lead 
to a rapid elimination of the residuals and require little initiative 
on the part of the computer familiar with the sequence of step-by-step 
operations. The method, however, is not efficient in the case of 
panels with many bays in the direction of the stringer loads or panels 
with sheet covering of large shearing rigidity, since large forces are 
then Introduced Into adjacent stringers when one stringer is balanced. 
These forces in turn must be liquidated in successive operations with 
the consequence that the procedure becomes time consuming. Also in 
problems involving many loading conditions it may be expeditious to use 
the electric -analogy method described in the section entitled 'the 
Electric Analogue, " since in the relaxation method each new loading 
requires new step-by-step operations. 

In this section panels are discussed which are not excluded from 
application of the relaxation method by the foregoing considerations. 
They may be classified according to the boundary conditions of the 
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stringers into four groups.. Recommendations for each, group follow with 
a fifth subsection added containing suggestions for panels in which 
transverse forces and displacements are considered. 

fa) Panels with boundary conditions at both ends of stringers speci- 
fied in terms of force . - The following two procedures are reco mme nded for 
liquidating the residuals on a panel of this group: 

First procedure: 

1. Consider each stringer isolated by cutting the sheet and the 
transverse reinforcing elements. Select the stringer for which the 
algebraic sum of the external forces is the largest. Displace the entire 
stringer as a rigid body (block displacement) until this sum vanishes. 

2. Balance one end Joint of the stringer by displacing the adjacent 
Joint on the same stringer. 

3* After step 2 is completed the end Joint is balanced but the Joint 
that was moved is unbalanced. Displace the third Joint on the same 
stringer until the second Joint is balanced. 

4. Continue the procedure until the second end Joint is moved. In 
this last step both the end Joint and the adjacent one will be approxi- 
mately balanced at the same time since the algebraic sum of all the forces 
acting upon the stringer was zero after completion of step 1 and this 
equilibrium has been disturbed only slightly by the shear forces trans- 
mitted by the sheet during the individual operations. 

5 . Stringer 1 is now approximately balanced. Carry out the same 
procedure with the other stringers of the panel successively. 

6 . When all the stringers are approximately balanced, return to the 
first- stringer and balance it again by undertaking steps 1 to 4- Repeat 
the procedure with the other stringers -until all the residual forces can 
be considered negligible for engineering purposes. 

Second procedure: 

1. Consider each stringer isolated by cutting the sheet and the 
transverse reinforcing elements. Select the stringer for which the 
algebraic sum of the external forces is the largest. Displace the entire 
stringer as a rigid body (block displacement) until this sum vanishes. 

2. Displace one end point of this stringer so as to balance the 
residual thereon. 

3 . Displace by equal amounts the adjacent- Joint on the same stringer 
and the end Joint which was balanced in step 2 so as to balance this 
second Joint. The equilibrium of the end Joint will be disturbed only 
by a small amount due to shear in the sheet. 
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ij-. Displace "by equal amounts the. third Joint on the Bame stringer 
and the two Joints that were placed In approximate balance by the oper- 
ation described in step 3 so as to balance this third Joint. 

5 . Continue this procedure until the Joint next to the mid Joint of 
the stringer is balanced by equal displacements of all the Joints 
situated between it and the end Joint first displaced. 

6 . Eepeat the process described in steps 2 to 5 , starting from the 
other end Joint of the stringer and continuing to the midjoint from this 
direction. After this step is completed this stringer will be in 
approximate balance, the only residuals being those introduced by shear 
in the sheet. 

7 . Consider next the stringer on either side of the approximately 
balanced stringer. Undertake a block displacement so as to equilibrate 
externally the stringer under its residual forces. 

8 . Start at one end Joint of this stringer and apply steps 2 to 6 . 
This second stringer will be placed in approximate balance thereby, 
while the balance of the first stringer will be disturbed only through 
the shear in the sheet. 

9 . Either return to the first balanced stringer or proceed to the 
next stringer on the other side. Each newly considered stringer is 
first externally equilibrated under the external and residual forces by 
a block displacement. Then from each free end the residuals are 
balanced by group displacements involving equal displacements of all the 
Joints situated between the one in question and the free end. Continue 
to balance Individual stringers until a.ll are balanced. 

The relaxation tables for the panel shown in figure 1, for which 
table 1 is the operations table, are used to demonstrate the first 
and second procedures and . are given as tables 2 and 3 , respectively. 

It will be noticed that this operations table considers the displacements 
of only the Joints on the left half of the panel. The panel is symmet- 
rical and is symmetrically loaded. Therefore, the displacements in the 
balancing process are undertaken symmetrically and only those of the left 
side Joints need be considered, those of the right being correspondingly 
equal. Since this panel has only three bays along each axially loaded 
stringer, the internal balancing process is undertaken from one end of 
the stringer only. 

ft ) Panels with boundary conditions at one end of stringers specified 
in terms of force and at other in terms of displacement . - This type of 
problem occurs, for instance, when one end of the panel is attached to a 
rigid body which is either held fixed in its position or is displaced a 
given amount. The recommended procedure for panels of this group is the 
same as the second procedure for panels in case (a) with two exceptions: 
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(l) Wo "block displacements are needed (or possible) to equilibrate the 
stringers externally and (2) the internal balancing process can be 
started only from the one free end of each stringer. 

The method is demonstrated on the panel shown In figure 2 . It is 
identical with the panel used for case (a) with the exception of the 
fixed lower ends of -the vertical stringers. The operations table is 
identical with that— of the previous panel except that no block and 
no v^ and v Q displacements are admissible. The relaxation table is 

given as table 4 * 

(c) Panels with boundary conditions at~both ends of stringer 
specified in terms of displacement . - Experience on panels of this type 
indicates that— although no systematic process of balancing the residuals 
can be recommended, the direct relaxation process is rapidly convergent. 

By starting from the midpoint Joints on a stringer and by balancing 
successive Joints toward the two fixed ends, the equilibrium position 
can be approximated rapidly. A further suggestion regarding this type of 
panel is contained in the later section 'If ilea Tables. " 

(d) Panels with irregularly specified boundary conditions . - For such 
panels a combination of the methods discussed under cases (a), (b), and 
(c) is rec ommen ded . By Judicious use of block and group operations 
similar ter those of cases (a) and (b) rapid convergence of the relaxation 
procedure will ‘be obtained. 

(e) Panels in which transverse displacements and forces are 
considered . - There are two general procedures for treating panels in 
which the transverse displacements and forces, visually considered negli- 
gible, are treated. These are described in the following paragraphs: 

First procedure: 

The procedure discussed under cases (a) and (b) can be applied to 
panels with cut-outs. The stringers are approximately balanced in the 
direction of the major axial forces by these procedures and then the 
residuals noraal to this direction are considered. The same step-by-step 
operations can be applied in balancing transverse stiffeners under these 
transverse axial forces. The process of first balancing the stringers 
in one direction, then balanc in g the stiffeners in the normal direction, 
and then returning to the originally balanced stringers will be quite 
rapidly convergent for panels with sheet of low shearing rigidity. 

Second procedure : 

For panels with cut-outs requiring consideration of the transverse 
forces another procedure, which is demons trated in reference 4 , can be 
used. The panel Is first considered to have continuous sheet and stringers, 
as if the cut-out did not exist, and the displacements for equilibrium of 
this panel under the external loads are determined by the usual methods. 
These displacements are then applied as a first approximation to the 
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distorted shape of the actual panel with cut-outs. Displacements leading 
to a closer approximation are then undertaken. This procedure is found 
to he reasonably successful for the cases investigated in reference 


Matrix Calculus Method 

The operations table together vith the external forces can be 
considered as a system, of linear equilibrium equations vith the magnitudes 
of the displacements as the unknowns. Therefore, the methods of matrix 
calculus can be applied to find the solution of this system by direct 
mathematical means. The method described in reference 9 is reco mme nded 
since a check on the calculations is maintained at each step in the 
process of solution. 

Matrix methods of solution have several advantages. After the 
operations table is established by trained engineering personnel, the 
solution can be obtained by computing personnel familiar vith the matrix 
calculus method. Under some conditions this economic advantage may be 
4 important. For reinforced panels vith sheet of high shearing rigidity 

the relaxation procedures are slovly convergent even vhen the reconmen- 
dations given in the preceding section are observed. The matrix calculus 
* method is not affected by this physical characteristic of the structure. 

When the number of equations is greater than 30 or 40, the vork of 
computation becomes inconveniently large. Therefore, for panels having 
a sheet covering of small shearing rigidity relaxation methods are 
recommended. When the sheet covering is very rigid in shear the matrix 
method is likely to be more advantageous because the routine operations 
of the matrix method can always be carried out if enough time is allowed. 

The equations of equilibrium for the panel shown in figure 1 are 
given by table 1 and are presented as follows to illustrate how the 
operations table and the external forces can be considered as a system of 
equilibrium equations: 



- 55-2t b + 2-00v e + 51.2r p = 0 

2-00v B - 101-6 v e + 4.00v p + 46.8 tj + 2.00Vg. = 0 

51.2vj + i»-.00v E " + 2.00Vj + 51.2v.g- = 0 

46.8v e + 2.00v p - 1D1.6 vj + ll-OOVg. + lf6.8vjj + 2.00v 0 = 0 

2-OOv.g + 51.2v p + If.OOVj, - 110. Wg + 2*00 Vjj + 5l.2v Q = 0 


46-8vj. + 2.00Vg - + 2.00v 0 + 60 x 10^ = 0 

2.(X)Vj + 51.2Vg. + 2.00Vjj - 55*2Tq + 60 X 10^ = 0 


5 




(1) 


-J 


In considering the operations table and the external forces as a system of eqidllbrlum 
equations, care must he taken to restrain, enough Joints so that the position of the panel as 
a rigid body is fixed. In the present case v A and Vjj are assumed to he zero, and since 

only displacements in the y -direction are considered in this problem, this restraint is 
sufficient. 
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Growing -Unit ' Method 

For reinforced panels with sheet of high shearing rigidity or with 
a large number of "bays in the direction of the axial forces, the relaxation 
procedure is not rapidly convergent. In such problems either the matrix 
calculus or the growing -unit method is rec omme nded. The latter can "be 
applied only to panels the boundary conditions of which are specified in 
teims of force at least at one end of the stringers. 

The growing -unit method applied to reinforced panels is as follows. 

The joint at the free end of an arbitrarily selected unbalanced stringer, 
called hereinafter the principal joint and the principal stringer, 
respectively, is displaced so as to liquidate the residual on this Joint. 

At the same time the Joints lying on adjacent parallel stringers and the 
same transverse stiffener are displaced so that the residuals that would 
be otherwise introduced by shear from the balancing of the principal 
Joint as well as any external forces applied to these joints are likewise 
liquidated. In the second operation the next joint on the principal 
stringer is relaxed while the previously balanced joints on the first 
transverse stiffener and the Joints on the second transverse stiffener 
are kept in balance by suitable displacements- After this second 
operation no residuals remain at the joints of the first two transverse 
stiffeners . After a sufficient number of repetitions of the procedure 
all residuals will be confined to reaction points or will be liquidated; 
the panel will then be in equilibrium. 

This procedure is demonstrated on the panel shown in figure 3 • The 
physical properties of the panel are the same as those of the previously 
discussed panels except for the additional bay in the direction of the 
axial forces. Actually the convergence of the relaxation method for this 
panel would be quite rapid, but for convenience the growing-unit method, 
applicable when this convergence is slow, is demonstrated thereon. 

Table 5 Is the operations table for this panel and contains not only the 
individual operations but also the group operations of the growing -unit 
method- Table 6 Is the relaxation table in which these group operations 
are used. 

The group operations given in table 5 require some explanation. In 
order to avoid Introducing a Y^ residual when Joint A is relaxed by 

application of operation ( 1 ), a v-g displacement is applied, the magni- 
tude of which can be calculated from the equation 


- 55-2Vg + 2.00 =0 (2) 

Thus operation ( 9 ) is v.g = ( 2 / 55 . 2 ) = O.O 362 and (10) Is a group 

operation equal to the sum of operations (l) and ( 9 ), which liquidates 
the residual Y^ without introducing a Y-g unbalance. 
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After operation (10 ) Is used, unbalances erist at- joints E and F, 
that is, on the second transverse stiffener. In order to balance these 
without disturbing the recently established balance at A and B, two 
group operations are developed: one permitting the balancing of~E and 

one permitting the balancing of F. The magnitudes of v A and v-g 

required to maintain the balance of A and B when a displacement 
of v-g - 1 is undertaken are given by the following equations: 

- 50. 8v A + 2-oOv-g + 46.8 = 0 


(3) 


2.00v a - 55>2v-g + 2.00 = 0 


These are satisfied by v A = 0-921, operation (11), and v-g = 0 . 0695 , 

operation (12). Operation (13) is therefore established as the sum of * 

operations (3)> (11), and (12) . The magnitudes of v A and v-g 

required to maintain the balance of A and B when a displacement - 
of v_ = 1 is undertaken are given by the following equations: 


- 50 .8v a + 2 . OOv-g + 2.00 = 0 




y 


W 


2.00v a - 55.2v b + 51.2 = 0 


These are satisfied by v A = 0.0758, operation (14), and Vg = 0.923, 

operation ( 15 ) . Operation (l 6 ) is the sum of operations (4), (l4), 
and ( 15 ). Since group operations ( 13 ) and (16) both introduce y e 

and Y-g, forces, the magnitudes x^ and x^ of these groups 

required to liquidate the -111 -pound and - 9 -pound residuals at E and F, 
respectively, are given by the following equations: 


- 58-3x 13 + 9-^ l6 - ill * 0 

> 


9.4 x ^ - 62. - 


(5) 


9 =» 0 
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Thus = - 1*975 and = - 0.444. Joints E and F are "balanced 

■without disturbing the balance of A and B by the use of these multiples 
of operations ( 13 ) and (l 6 ). 

In eliminating the residuals at joints J and E multiples of 
operations ( 13 ) and (l 6 ) are applied since these operations permit 
displacements of E and F to be undertaken while the balance at A and B 
is left undisturbed. When Joint J is displaced a unit amount, multiples 
of operations ( 13 ) and (l 6 ), defined by the following equations, are 
used so that the balance at A, B, E, and F is maintained: 

\ 

► (6) 


58 . 3 X 13 + 9 *^ i6 + 46.8 = 0 


9>4x - 62 . 6 x ^ + 2.00 = 0 

13 Id 


The solution to these equations is x-^ = O. 828 , operation ( 17 ), and 

= O.I58, operation (l8). Operation ( 19 ) is the sum of operations (5), 
( 17 ), and ( 18 ). 

In a similar manner all the individual and group displacements 
described in table 4 are found. It may be mentioned that in the present 
example no shearing stresses were set up in the middle bays because of 
the symmetry of structure and loading. The original operations table 
was already established in a manner which complied with these require- 
ments of symmetry. When such Is not the case or when there Is a greater 
number of stringers in the panel, displacements of principal stringer 
Joints will, in general, cause residuals to appear at more joints so 
that three or more, rather than two, simultaneous equations have to be 
solved at each step. 


Niles Tables 

In reference 10, A. S. Niles demonstrates for the solution of rein- 
forced-panel problems a method which essentially parallels the previously 
described relaxation method. The Niles method Is a procedure for 
balancing a stringer by the use of tables which give the displacements 
of each joint on the stringer required to liquidate a residual on a given 
Joint of the stringer. The tables are worked out for various end 
conditions and sheet shearing rigidities. 

Since reference 10 contains tables only for sheet of relatively low 
shearing rigidity, the Niles method is limited in this respect in the 
same way as the relaxation method. However, the tables can be employed 
on stringers with the boundary conditions at both ends specified in 
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terms of displacement} for such problems no step-by-step routine 
relaxation method has been recommended. Also by use of the tables exact 
balance of a stringer is gained after a single displacement of each 
Joint, whereas in the relaxation method, because of the shear, small 
unbalances remain after each Joint is moved. 

On the other hand, the relaxation method can be applied to stringers 
with irregularly spaced Joints for which no tables were set up by Niles. 

Since in reference .10 several examples of the procedure-are given, 
no application of the Niles method is shown herein. 


Electric Analogue 

Another convenient method of solving the problem of force distri- 
bution in a reinforced panel is that^in which the voltages are measured 

in an electric network which is so constructed as to make it a complete 

analogue of the reinforced panel. "When suitable electric equipment is 

available, an analogous network can be hooked up and tested with very 

little work- A particularly attractive property of - the stress -analysis 
procedure by means of electric measurement is the ease with which the 
effect upon the stress distribution of changes in loading and in dimen- 
sions of the various structural elements of the reinforced panel can be 
investigated. This permits the development of an efficient=deeign with 
little analytic work. 

The analogy between the forces transmitted through the different- 
structural elements of the reinforced panel and the currents flowing 
through the various branches of the direct -current network can be 
explained with the aid of figures 4 and 5. The problem investigated is 
the so-called "one -dimensional shear lag." It is assumed that the trans- 
verse stiffeners are infinitely rigid so that the vertical, or longitu- 
dinal, displacements v alone need to be determined. The portion of 
the sheet - covering considered effective in tension or compression is 
added to the cross-sectional area of each stringer and the panels of 
sheet are assumed to carry shear stresses only. A consequence of these 
assumptions is that the shearing stress must be constant in each panel. 

The analogous direct -current network contains as many binding posts 
as the number of Joints in the reinforced panel. Adjacent - binding posts 
are connected by conductors having prescribed resistances R. Prede- 
termined electric currents I, which correspond to the forces F 
applied to Joints A and B of the reinforced panel, are introduced into 
the network at points A and B. 

Xt is now recalled that in the relaxation method the joints of- the 
panel are first assumed to- be rigidly fixed to a rigid wall behind the 
t panel. The external loads are first applied, to these .rigid pegs, 
referred to as the "constraints. " The panel is obviously in equilibrium 
under these conditions but this artificial equilibrium is entirely 
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different from that prevailing in the actual panel, -which is not attached 
to any rigid vail. The actual state of equilibrium is approached by the 
step-by-step procedure of the relaxation method, in each step of which 
one single constraint is removed and the corresponding Joint is displaced 
until it reaches its equilibrium position in the system in which all the 
other Joints are still rigidly fixed. 

For instance when Joint 1 of the reinforced panel is moved through 
a distance v in the positive direction, this displacement imposes 
forces upon all the adjacent Joints numbered from 2 to 9. Three typical 
forces are given by the equations: 



( 7 ) 


(8) 


( 9 ) 


where 

F 8l , F gl , Fg^ the forces acting upon Joints 8, 9, and 6, respectively, 

because of the displacement of Joint 1 

E modulus of elasticity of stringer 

G- shear modulus of sheet 

t thickness of sheet 

v displacement of Joint 1 

In the case of the analogous network it can be assumed that the 
potential of each binding post is zero at the outset- If there Is no 
potential difference, no current flows between the posts. It can be 
imagined that the currents introduced at points A and B are taken out of 
the system by means of some imaginary conductors • However, the actual 
distribution of currents in the network prevails without the aid of the 
ima ginary conductors. This actual state can be approached also by means 
of a step-by-step, approximation -type calculation. For instance it can 
be assumed first that the potential of binding post 1 is elevated to the 
value V. After this change there is a potential difference between 
binding posts 1 and 8 and consequently a current will flow from post 1 
to post 8. The magnitude of this current can be calculated from the 
equation 


*81 “ V E 8l - °8l T 


(10) 
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where Rg^ is the resistance and Cg^ = -y ®8i is tile conductance of the 

conductor "between posts 1 and 8. Similarly tlie current flowing from 
postrl to post 9 is 



= C 9 1 V 


(ll) 


The current flowing from post 1 to^post 6 is 



= C 61 T 


( 12 ) 


Comparison of equations (7) to (9) with equations (10) to (12) 
reveals an analogy between the effects of a displacement v of Joint 1 
and the raising of the voltage of "binding post 1 "by an amount V. The 
current- caused "by the change in potential corresponds to the force 
caused "by the displacement, provided that the conductance of each 
conductor is made equal to the influence coefficient in the corresponding 
force equation. Hence 


°8i = 

EA Gbt 
b " 2a 

( 13 ) 

V 

= qbt- - 
4 a 

( 14 ) 

c 6i 

_ Gbt 
2a 

( 15 ) 


In the relaxation procedure the equilibrium state is approached by 
displacing individually the Joints and summing the effects of each 
displacement. In exactly the same way the actual distribution of the 
currents in the network can be determined by changing individually the 
voltages of each binding post and summing the effects of these changes. 
In the reinforced panel equilibrium, is obtained when at each Joint the 
sum of the external forces and of all the internal forces caused by the 
displacements Is zero. The forces are considered positive If they are 
directed as the positive displacements. In the form of an equation. 
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An. analogous equation in the direct -current network is furnished by 
Kirchhoff 's first law, according to which the sum of the currents flowing 
into any "binding post must "be zero- Currents in the direction of any 
"binding post are considered as positive- In the form of an equation, 

]Ti = 0 (17) 


Comparison of the last two equations reveals that the conditions of 
equilibrium for the reinforced panel and EIrchhoff 's first law in the 
case of the direct -current network complete the analogy of the two 
systems considered. It Is possible therefore to construct an electric 
network with the same configuration of binding posts as that of the Joints 
of the re inf orced panel. The conductances of the conductors connecting 
the binding posts must be so chosen as to make them proportional to the 
corresponding influence coefficients in the operations table of the 
reinforced panel- If then currents are introduced at the binding posts 
which correspond to the Joints at which external loads are applied, the 
distribution of the currents in the network will be the same as the 
distribution of the forces between the various structural elements of 
the reinforced panel- 

In the first applications of the relaxation process to reinforced 
panels each Joint was displaced until equilibrium was established. It 
was noted in the section dealing with the solution of the problem by 
matrix methods that this procedure permitted rigid body displacements of 
the structure. Blgid body displacements can be eliminated If one or more 
Joints are considered as rigidly fixed- In the case of the reinforced 
panel of figure k the degree of freedom of motion of each Joint is one, 
because the problem is considered as a one -dimensional shear lag problem. 
Consequently It suffices to fix one single Joint so that it is prevented 
from displacing vertically. However, If Joint C, for instance, is fixed, 
the symmetry of the structure and loading requires the simultaneous 
fixation of Joint D. 

In the analogous network binding posts C and D are given predeter- 
mined values of the potentials by connecting them to the ground. It is 
customary to attribute the value zero to the potential of the ground- 
Consequently V c and Y D are zero Just as in the reinforced panel 

v c and v-p are zero- 

It will be noticed that in figure 4 the direction of F at 
Joints A and B is upward, whereas the direction of I at binding posts A 
and B in figure 5 is downward. This corresponds to the difference in 
the sign convention in the two systems- In the panel upward forces were 
considered positive and in the network currents flowing toward the 
binding posts were given the positive sign- The directions of the forces 
and the currents at points C and D are the same- This again corresponds 
to the correct signs required by the sign convention since the downward 
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forces at these points are negative Just as the currents which flow away 
from the binding posts are negative. Hence the reinforced panel is 
under the action of external tensile forces, whereas through the network 
currents are flowing in the downward direction. 

In the case under discussion it is easy enough to introduce the two 
equal currents at posts A and B and to- regulate their magnitude by means 
of an adjustable rheostat. However, when there are a number of impressed 
currents of different magnitude stipulated, their adjustment may become 
a lengthy trial -and -error procedure. In such cases it is advantageous 
to employ a number of commercially available electronic devices, known 
as constant -current generators, which have the property of maintaining 
a constant current independently of the properties of the network. 

When the construction of the network Is completed and the required 
external currents are introduced, the deflection of any Joint of the 
reinforced panel can be obtained by measuring the potential of the 
corresponding post in the network with respect to the ground. This 
quantity multiplied by the conversion factor Is the relative displacement 
of the corresponding Joint of the reinforced panel with respect" to - the 
fixed points C and B. In most cases, however, the displacement quantities 
are of interest only Indirectly and the main quantities sought are the 
forces In the stringers and the shear stresses in the sheet. These 
quantities can be obtained in a simple manner by multiplying potential 
differences by the appropriate conductances and by the conversion factor. 

For instance when the force in stringer segment 1-8 is sought, the 
voltage drop between posts 1 and 8 must be measured and multiplied by 
the conductance Cg^ and the conversion factor. This is a consequence 

of equations (7) and (iO ) . Similarly when the shear stress in panel 1689 
is required, the voltage drops in conductors 1-6 and 8-9 have to be 
measured. From figure 4 the average displacement of stringer segment 6-9 

is + Y g)j^- and the average displacement of stringer segment 1-8 

Is ( v 3_ + vg)/ 2 . The difference of these two average displacements 

multiplied by Gtb/a Is the shear force transmitted from, the panel 
to stringer segment 6-9. Consequently the sum of the displacement differ- 
ences Vg - v^ and v^ - Vg multiplied by the influence coefficient 1-6 

is the shear force sought- In other words the sum of the voltage drops 
from post 1 to-post 6 and frcaa post~8 to post 9 multiplied by the 
conductance Cgj_ and the conversion factor is the shear force in question. 

This shear force divided by the length b gives the average shear flow in 
panel 1689 and this shear flow divided by the thickness of the sheet is 
the average shear stress. 

With the cooperation of the Department of Electrical Engineering a 
network was constructed at the Polytechnic Institute of Brooklyn which 
was the analogue of the reinforced panel investigated earlier at P 3 BAL 
both experimentally and by relaxation methods- The results of these 



NACA TN No. 1786 


19 


investigations are described in reference 3 . The constant currents were 
introduced by means of constant -current generators. In the electrical 
system the unit of the potential was chosen as 1 volt and that pf the 
current as 100 mllliamperes . Then the unit of the conductance had to he 
a millimho and that of the resistance, a kilohm. In the mechanical 

system the unit displacement was 10 "^ inch and the unit force, 1 pound. 
Consequently in this problem the voltage differences had to he multiplied 

hy the conversion factor 10 inch per volt in order to obtain displace- 
ments- The factor converting currents into forces was 10 pounds per 
ampere. The results of the measurements were in excellent agreement with 
the results quoted in reference 3 . 

Similar experiments were carried out by R. E. Newton and M. E. Engle 
at the Curtiss -Wright Corporation, Airplane Division, in St. Louis and 
are described in two reports listed as references 11 and 12. Newton's 
approach to the problem is fundamentally the same as the argument given 
herein. However, his electric network is slightly simpler since it does 
not contain the conductors arranged diagonally in the system shown in 
figure 5 . 'The network of figure 5 was chosen in this report in prefer- 
ence to Newton 's simpler network since by this presentation the iden- 
tity of the conductances of the network and the influence coefficients 
used in the other parts of this report could be established. 

It should be mentioned that in many cases it is possible to construct 
a dual type of network in which the currents correspond to the displace- 
ments of the Joints of the reinforced panel and the potential differences 
correspond to the forces in the stringers and in the sheet covering of 
the panels. In this type of network the external loads can be introduced 
more easily as impressed potential differences- However, the network 
described herein Is more advantageous since it can always be constructed 
directly from the geometry of the reinforced panel. 

The usefulness of the analogue with the direct -current network 
breaks down when the influence coefficient in equation ( 7 ) becomes 
negative. In such a case the conductance and consequently the resistance 
of the corresponding branch of the network should be negative] this is 
obviously impossible. However, the situation can be usually remedied in 
the case of one -dimensional shear lag problems. The fundamental assump- 
tions of the problem are not changed if a number of additional horizontal 
bracing elements are introduced in the panel since all of them are assumed 
to be infinitely rigid. If, however, the panel length b Is reduced 
to one -half its original value, then the negative term in the influence 
coefficient appearing in equation ( 7 ) Is halved and the positive term is 
doubled. In most cases this will suffice to change the sign of the 
influence coefficient. When such Is not the case distance b can be 
reduced in any other suitable ratio. 

Negative influence coefficients can be realized If the analogous 
network is fed by an alternating current. The quantity corresponding in 
an alternating -current circuit to the resistance of the direct -current 
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circuit Is the Impedance • In the Impedance the inductance retards the 
phase of the current and the capacitance advances it so that" the two 
have opposite effects. If one Is designated as positive, the other Is 
negative- However, no Inductance is entirely free of resistance and for 
this reason the accuracy of a complicated alternating -current-network 
may not' he sufficient for the solution of some of the problems encoun- 
tered In practice- 

The use of the electric analogue for solution of shear lag problems 
is recommended when several similar panels with many loading conditions 
are to be analyzed- For such a case the construction of the analogous 
network, the variation of the loading by varying the impressed currents, 
and the determination of the potentials at the binding posts would be 
simpler than any analytic method of solution. — 


FOSELAGE KINGS 
Introduction 


In reference 5 numerical procedures for the determination of the 
bending -moment— distribution in fuselage rings are developed and demon- 
strated on several simple and internally braced fuselage rings- The 
number of redundant internal bracing elements increases little the work 
involved in establishing the operations table for the ring and affects^ 
not at all the amount of numerical work in the solution of- the operations 
table- This nonsensitivity to the number of redundances constitutes the 
advantage of this method in the analysis of fuselage rings. 

The methods suggested for the solution of- the system of equations 
represented by the operations table and the external forces are three: 
relaxation, matrix calculus, and growing -unit. The latter two may be 
considered as direct mathematical methods and as in reinforced-panel 
problems require only computing personnel. For the analysis of isolated 
fuselage rings of coup lex shape the use of these direct methods is 
recommended since an accurate solution is assured in' a reasonable length 
of time, whereas the relaxation method may notrlead to"sufficientiy 
accurate results even after considerable effort has been expended- How- 
ever, for simply shaped rings and for problems of stress distribution in 
sheet-- stringer, and ring combinations, application of _the relaxation 
method to fuselage rings is advantageous . For this reason the relaxation 
method for fuselage -ring problems is presented and new, more rapidly 
convergent procedures are developed. 

It has notrbeen found possible to maker concrete reommmendations for 
relaxation procedures which are rapidly convergent for all types of ring 
and loading. However, satisfactory procedures for several distinct types 
of ring and loading are demonstrated and explained in some detail. It is 
felt that consideration of these examples will suggest to the analyst 
means of solving more rapidly other ring and cylinder problems which are 
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not efficiently attacked "by direct mathematical means. The procedures, 
which involve essentially appropriate combined operations, are demon- 
strated on two rings solved in reference 5 "by the usual relaxation 
methods and on a new internally braced ring. Application of the growing-: 
unit and matrix calculus methods to the latter problem is made to demon- 
strate these methods and to verify the results of the relaxation 
procedure . 


Torsion of a Circular El ng 

In reference 5 the bending -moment distribution for a simple 
circular ring with antisymmetric loading consisting of concentrated 
forces and distributed and constant shear flow is determined by appli- 
cation of n um erical methods . The di-ma ns ions and loading for this ring 
are shown in figure 6 and the operations table Is given as table 7. 
Belaxation methods are applied to the solution of this ring problem in 
reference 5. By a process of Increasing si 1 the residuals in such a 
proportion that one key operation would liquidate them, all to within 
the desired degree of accuracy, the residuals were reduced to within 
2 percent of the maximum applied load in 12 operations. 

In the present report combined operations which increase the 
rate of convergence are demonstrated. Tangential and angular dis- 
placements of A and C balance these points in four operations and place 
all remaining residuals at B . Since no tangential forces exist at 
A and C, the force residual at B must be vertical and the moment 
residual, equal to the couple of the vertical forces. Suppose the 
residual moment at B is liquidated by a rotation of that Joint while 
the balance of A and G is preserved by suitable displacements of A and C. 
Then from equilibrium considerations the residual forces at B must also 
be liquidated. Thus in five operations balance will be obtained. This 
procedure is used and proves to be satisfactory . 

In order to b alan ce the residuals at A two combined operations 
are developed. The first combines a unit angular displacement w^ 

with a tangential displacement u^ such that no tangential force at A 

results when the two individual operations are simultaneously applied. 

The forces and moments introduced by the individual operations as well 
as by the combination are given in the following table: 



Forces and 
moments 
Operation ' — 

*A 

t a 


— 

— 

*c 

T c 

v, = 10 "3 radian 
A 

-281.95 

-49.079 

-29.966 

- 4-733 

64.675 

0 

0 

u A - - 0 . 9384 fi x 10 " 3 In. 

46.060 

49.079 

-60 . 696 

21.060 

-48.347 

0 

0 

— » Operation A = 1 

- 235-89 

0 

-90.662 

16.327 

16.328 

0 



0 


The second operation combines a unit tangential displacement u A with an angular rotation v A 

such that at A no moment arises from the combined operation- The forces and moments introduced 
by the individual operations as -Hell as by the combined operation are given in the following 
table: 


Forces and 
momenta 

Operation 

*A 

t a 

% 

% 


*0 

T c 

u^ = 10 ‘3 in. 

49.079 

-52.296 

64.675 

- 22.441 

51.516 

0 

0 

w A * -0.17407 x 10 “ 3 radian 

49.079 

8.5432 

5.2162 

0.82387 

-II.258 

0 

0 

y~ -> Operation B ■ 1 

0 

43-753 

69*891 

-21.617 

40.258 

0 

0 



























Thus "by using the necessary amounts of the combined operations A and B Joint A is balanced 
in two operations. Two similar operations are found for Joint C and are given, as follows with- 
out explanation: 


Forces and 
, '-~^moments 
Operation 

*A 

t a 




K c 

- ‘ 

T 

c 

v c = 10 "3 radian 

0 

0 

56.512 

8.842 

6.632 

- 157.899 

-1.563 

Up = - 4 . 85^0 x 10 “3 1^. 

0 

0 

- 32.192 

-e.5435 

-0.33250 

7.5668 

1-563 

y~ —> Operation C = 1 

0 

0 

24.320 

6.2985 

6.2995 

-150.31 

0 


Forces and 
^^^-■^momentB 
Operat ion 

"a 

t a 

% 


*6 

*C 

T o 

Up = 10 "3 in. 

0 

0 

6.632 

0.524 

0.0685 

- 1.563 

-0.322 

w c = -O.OO98987 x 10 "3 radian 

0 

0 

- 0-55939 

-0.08752 

-O.0656 

1-563 

0.01547 

-> Operation D = 1 

0 

0 

6.0726 

' 0.43648 

O.OO29 

0 

-0.30653 


In order to balance the residuals at B without disturbing the balance at A and C obtained by 
use of operations A to D, combined operations involving tangential and angular displacements of 
A and C and a unit rotation of B are developed. If Joint A is to remain in balance when a rotation 
of B is undertaken, Joint A must be rotated and displaced in such a manner that the tangential 
force and the moment introduced at A by this rotation of B are .equilibrated* Since the angular m 
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displacement introduces tangential forces at A, and the tangential 
d i s plac ement Introduces moments, two simultaneous equations must "be 
solved for the unknown tangential and angular displacements. The 
equations for A are: 


- 28i.95v a - k9-079u A - 2 9.966 x 10 "3 = 0 

- 49«079w a - 52 . 296u a + 64.675 x 10 " 3 = 0 


> 


(18) 


The solution to these equations is V A = “ 0-38434- X 10 “8 radian and 

u A = 1.5974 X 10 “3 inch. A unit rotation of B and tangential and 

angular displacements of C are combined in equations ( 19 ) so that the 
tangential force and moment introduced at C by the comb ine d operations 
are zero. 


- 157 • 899v c - 1.5631^ + 56.5117 x IO ' 3 = 0 
- 1.563 v c - 0-322 uq + 6.632 x 10 "3 = 0 


(19) 


The solution to these equations is vq = O.I6180 X 10 "3 radian and 
Uq = 19.811 x 10 “3 inch. 


If the forces and moments introduced by the three sets of 

displacements (unit rotation of B, the tangential and angular dis- 
placements of A, and the tangential and angular displacements of C) 
are combined, a combined operation is obtained such that only forces 
and moments at B and radial forces at A and C are introduced. These 
latter forces are of no interest in the relaxation procedure since 
they are equilibrated automatically by the other half of the ring. The 
combined operation from these three sets of displacements is given in 
the following table: 



" — Forces and 
^^^-^moments 
Operation 

*A 

t a 

% 

% 

% 

ffc 

T c 

vtg = ID "3 radian 

-29.966 

64.675 

-439.849 

31.443 

-50.642 

56.5117 

6.632 

V A = x 10 “3 radian 

IO8.37 

18.863 

11.517 

1.8191 

-24.857 

0 

0 

u A = 1.5974 X 10-3 in. 

- 78.399 

-83.538 

103.31 

-35.047 

82.292 

0 

0 

vq = O.I6180 x 10 “3 radian 

0 

0 

9.1436 

1.4306 

1.0730 

-25.540 

-0.25289 

Uc = 19.811 x 10 “3 in. 

0 

0 

131-39 

10.381 

1.3570 

- 30-964 

-6.3791 

y~ Operation E = 1 

0 

0 

-184.49 

9-2267 

I— 

9*2230 

0 

0 


The relaxation table using these five combined operations, A bo E, Is given as table 8. The 
balancing process was carried put on a slide rule and after five operations all the residuals were 
reduced to negligible quantities. From the magnitudes of these group operations the total individual 
displacements of A, B, and C can be found and the 'unknown, radial forces at A and C calculated. 

The procedure just described involves essentially the development of group operations so that 
full advantage of the symnetry properties of the ring may be realized. This method is applicable 
to other rings* The internally braced circular ring subjected to antisymmetric loads and analyzed 
in reference 5 can be treated in the same way as this simple r in g. If these rings had been sym- 
metrically loaded, the force residuals at B, after A and C had bean balanced by simple radial 
displacements, would have a horizontal resultant. By combining radial and tangential displacements 
of A and 0 such that the resultant force introduced at B is horizontal and such that A and C remain 
in balance, the horizontal resultant at B could be liquidated by application of such a combined 
operation- The moment residual at B 1 b not necessarily eliminated when the force residual at B is 
balanced. Joint B must be rotated while A and C are displaced radially so that the m o ment at B is ^ 
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liquidated and Joints A and C are kept in balance. If the process of 
liquidating first-the residual force and then the moment at B, preserving 
in each operation the balance at - A and C, is not rapidly convergent^ -- two 
equations for the equilibrium of B can be established and solved for the 
required amounts of the combined operations . 

Thus the foregoing procedures for both the symmetrical and antl- 
symmetrical loading can be applied to - any ring singly symmetrical with 
only one Joint between the center line of symmetry Joints. It may, 
therefore, be advantageous in some ring problems to combine Beveral bars, 
as in the method of the growing unit, such that only one Joint between 
the boundary Joints has independent degrees of freedom- This will permit- 
use of the foregoing procedure . 

Sufficient accuracy for most engineering purposes can be obtained 
in the computations of this procedure by the use of a slide rule through- 
out. Although the combined operations shown herein were obtained by 
the use of a computing machine carrying five significant figures, the 
procedure was first demonstrated with the use of- a slide rule for all 
calculations. The results of-the two sets of calculations are in good 
agreement, thus indicating the sufficiency of slide -rule accuracy. 


Egg-Shaped Ring 

Figure 7 shows the dimensions of, and loading on, a ring which is 
analyzed in reference 5 . The operations table for this ring is given 
as table 9* In this ring there are two points B and C between the center 
line of symmetry points A and D. By making the degrees oil freedom of 
either point B or C dependent on the other and on the adjacent center 
line of the symmetry point, one point with independent degrees of free- 
dom is established between A and D and the method discussed previously 
can be used. 

However, in order to demonstrate the simplicity and effectiveness 
of group operations, another approach is used. The center lines of 
symmetry points A and D are balanced by simple radial displacements of 
A and D. The midpoint of bar BC is assumed restrained tangentially so 
that only equal and opposite tangential displacements of B and C are 
undertaken. Because of the large extensional stiffness of bar BC as 
compared with the bending rigidity of the circular segments and because 
the ring is a lmo st symmetrical about a horizontal axis, such dis- 
placements of B and C liquidate approximately equal and opposite 
tangential residual forces at B and C, such as those which will be 
obtained at these points when the residuals associated with the other 
degrees of freedom are small. 

If the balance at A and D is preserved by appropriate combinations 
of the radial displacements of A and D with the required displacements 
of B and C and if the tangential residuals at B and £ are not— cons idered 
until the foregoing operation will li qui date them both, main attention 
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is focused on the radial force and moment residuals at B and C. In 
order to Balance these, no specific method of convergence is used hut 
the state of the residuals after each step is considered Before the 
next operation is selected. In this problem of egg-shaped rings and 
many other rings and in the complete cylinder problems this approach, 
utilizing physical properties of the system and eliminating or reducing 
extraneous forces and moments at each step in the relaxation process, 
may Be the most satisfactory method of solution. 

Table 10 is the relaxation table for the ring in question- The 
first tvo operations involve only radial displacements -which balance 
the 500 -pound forces at A and D. The largest residual then Is the 
radial force of 451 pounds at C. If point C Is displaced radially so 
as to balance this residual, a large moment and a large radial force 
are introduced at B. In order to reduce these extraneous forces and 
moments and to keep Joints A and D balanced, radial displacements of 
A, B, and D and a rotation of B are combined as shown by the following 
operations : 


- 3* 34833 v A + 8.922161*3 - 2.696l4v B = 0 


"7 


8-922l6v A - 327 . 866vg 
- 2 . 696 l 4 v A + 11.4697-1*3 


+ 11. 4697 v B + 8.10267 x 


- 4. 0099173 + 0.66158 x 10 


- 12 . 2400 ^ - 1.11900 X 10 - 1 *- 



( 20 ) 


0 



The solution of this system of equations is: v A = - 0 . 26384 x 10 "^ Inch, 

Wg = O.O3279 X 10 radian, v-g = O.43618 x 10 “^ inch, and 
v B = - 0.9024 X lO 4 inch. 


The forces and moments introduced by each of the individual operations 
and by the combination are given in the following table: 



28 


NACA TIT No. 1786 


4 


^ Forces and 

— ■ — ^moments 
Operation 

r a 

“b 



t c = 10 in. 

0 

8.10267 

O.66158 

0 

v A = -O.26384 x 10 in. 

0.88343 

-2.35^0 

0.71136 

- 1.0468 

w B = O.O3279 x 10 4 radian 

0.29258 

-10.751 

0.37612 

-0.42963 

! v B = 0.43618 x 10 in. 

-1.1760 

5.0028 

- 1.7491 

1.4984 

1 V D = -0.90242 x 10 in. 

0 

0 

0 

0 

y~ — > Operation F = 1 

0 

0 

0 

0.0219 

Forces and 
~~ -—foments 

Operation 

■ H c 

BC 

T c 

b d 

v c = TO” 1 *- in. 

-2.95622 

-1.90205 

-0.88929 

-1.1190 

v A =.- 0.26384 X 10 in. 

0 

0 

0 

0 

w B = 0.03279 X 10 radian 

-2.0082 

0.26570 

0 

0 

v B = O.43618 x 10 in. 

- 3 - 53^2 

0.28857 

0 

0 

v D = -0.90242 x 10 in. 

6.6350 

1.00981 

0.93626 

1.1190 

Operation F = 1 

-1.8636 

. -0.33798 

0.04697 

0 


The use of combined operation F is desirable in balancing the radial 
residual force at- C, since it also reduces the moment residual at C and 
adjusts the tangential residuals at B and C in the desired manner. 

The residual considered after use of operation F is Rg = 402 pounds. 
In order to balance it by a displacement while the balance at A 

is preserved, a v A displacement must-be undertaken as well. If 

v-n = 10 ~ k inch, then y. = 9 • = - 0 . 80522 x inch. The 

-3.3W33 x io 4 





































































forces and moments Introduced "by these individual operations as veil as by the combinations are 
given in the following table: 


Forces and 
Operation 

*A 

% 

% 

% 


V 

T c 

®D 

v B = 10 in. 

- 2.69614 

11.4697 

-4.00991 

3.4352 

-8.10267 

O.66158 

0 

0 

v A = -0.80522 x 10“^ in. 

2.6961 

-7.1843 

2. 1710 

- 3 * 19 ^ 

0 

0 

0 

0 

^>~ Operation G = 1 

0 

4.2854 

-I.8389I 

0.2403 

-6.102 (f\ 

0.66158 

0 

0 


Consider the effect of eliminating the Eg residual by use of operation G- The moment 
residual Eg would also be reduced by roughly 1000 inch -pounds, the Tg residual would be 
brought in closer agreement with the Tq residual, an 1^. residual of about 30 percent of the 
previous Eg residual of 451 pounds would be introduced, and a large Eg residual would be 

introduced. The last two effects are undesirable. However, by use of operation 3 T again, the 

residual can be balanced without introducing a new Eg residual- The large Bg residual is 

not so easily balanced unless a new combination Involving joints A, B, and D,is evolved- 

Suppose, therefore, that a rotation of C and a radial displacement of D are combined so that 
a moment at C can be eliminated and so that the balance of D is preserved by use of the combination. 
The individual operations and the combinations are given in the following table: 




Foroes and 
" — moments 
Operation^ 

% 

% 

% 


*c 


T c 

%> 

w c = 10 -if radian 

0 

-61.242 

-8.10267 

0 

-288.367 

-£.95622 

- 5*24667 

- 7.3524 

Up = -5-9294 x 10 "^ in. 

0 

0 

0 

0 

43*595 

6.63500 

6.1517 

7.3524 

— > Operation H = 1 

0 

-61.242 

-8.1027 

0 

-244.772 

3.6788 

B3 

0 


If operations G- and E axe combined bo that the moment at C introduced by the combination 
is zero, the resulting forces and moments are given in the following table: 


Forces and 
^^^^inoments 
Operation 

e a 




*c 


T 

C 

*D 

IX (0) 

0 

4.2854 

-I.8389I 

0.2403 

-8,10267 

0.66158 

0 

0 

- 0.033103 X (H) 

0 

2.0273 

0.26822 

0 

8.10267 

-0.12178 

-O.O2996 

0 

-> Operation 1 = 1 

0 

6.3127 

- 1.5707 

0.2403 

0 

0.53980 

-O.O2996 

0 


Use of operation I results in liquidation of the Rg residual, in reduction in the 
Eg residual, in adjustment of the Tg and residuals toward the desired equality, and in 
introduction of an residual of 138 pounds- The latter can be balanced by the use of 
operation F, which will preserve the balance of A and D and will not affect the Hg and Eg 
residuals . 
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After this fifth operation the T-g and Tq residuals are approxi- 
mately equal and opposite as desired- Therefore, a group operation, 
involving equal and opposite tangential displacements of B and C and 
sufficient radial displacements of A and D so that the latter remain 
Balanced, is developed in the following table: 


" Forces and 

moments 

Operation^'~~--~~^_^^ 

r a 

% 

% 

t b 

N C 

Ug = 10 -14 - in- 

3.96771 

-13.1011 

3.1352 

-30.9566 

0 

Uq = -10 in- 

0 

0 . 

0 

-26.2058 

5.21667 

v A = I.I85O x 10 in. 

-3.9677 

10.573 

-3.1919 

1.7017 

0 

t d = O.83669 X 10 In- 

0 

0 

0 

0 

-6.1517 

Operation J = 1 

0 

-2.5281 

0.2103 

-52.161 

-0.9050 

" — Forces and 

77 moments 

Operation ~~ — ^ 

% 

T c 

e d 



Ug = 10 in. 

0 

26.2058 

0 



Uq = -10 in- 

0.88929 

27.0833 

1.0375 



v A = I.I85O x ID"! in. 

0 

0 

0 

m 


Vjj = O.83669 x 10 in. 

-0.93626 

-0 . 86807 

-1.0375 



y~ -> Operation J = 1 

0.01697 

52.121 

0 




Use of operation J liquidates the T-g and Tq residuals and 

affects little the B alan ce in the other degrees of freedom. The remai ning 
residuals are considered negligibly small, the moment of 309 inch -pounds 
Being approximately 3 percent of the maximum moment in the ring- As in 
the previous problem the Individual displacements can Be determined from 
the magnitudes of the group operations and thus the unknown mome nts and 
tangential forces at A and p calculated- 
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Although the calculations of the group operations shown herein hare 
"been carried out on a computing machine with fire significant figures 
maintained wherever possible, sufficient- accuracy for engineering 
purposes can "be obtained by the use of a slide rule. In developing this 
procedure a slide rule was used for all computations and the results 
agreed satisfactorily with those shown herein. 


Oval -Shaped King with Internal Bracing 

The ring shown in figure 8 is used as a third example of the new 
relaxation procedures. As a check on the results of this procedure the 
system of equations given by the operations table and external forces 
is also solved by the exact mathematical methods of matrix calculus and 
of the growing -unit method. In order that the charts and tables of 
reference 6 could be used in determining the i nf luence coefficients, 
the following physical characteristics of the elements of-the ring are 
assumed: 

Segments AB and EF: 


7 = = 500 

1 = m = 0 . 10 

A 

0 = 45 ° 

6 2 

EE = 10° lb -in. 

L = 18.85 in. 

Segments BC, CD, and DE: 

7 = = 500 

£ = 0.10 

3 a 


0 = 30 ° 



HA.CA TN Wo. 1786 


33 


EX = 10 6 lb -in- 2 


L = 18.85 in. 


Segment EG-: 


7 = &E. = koo 
I 

I = ££ = 0.10 

A 


P = 0 

H = 10 5 It -in. 2 


L = 16.97 In* 


Because of the symmetry atout a line through AGE only one -half of the 
ring need te considered. Joints A, G, and E are then, restrained from 
rotating or displacing tangentially and cannot te subjected to radial 
forces. The assumed positive directions of the displacements and of 
the forces and moments at each joint are shovn in figure 8.- From the 
foregoing assumptions, the influence coefficients and the operations 
tatle given in tatle 11 are determined. 

The horizontal external forces of 1000 pounds at C and D are 
resolved into their tangential and radial components • Thus the 
external forces are: 


Kc = 9 ^ 5-93 lb 
T c = - 258.82 lb 

Rq = - 965-93 It 
T D = - 258.82 lb 


> 


( 21 ) 
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The matrix calculus solution of the system of equations given "by 
these external forces and by the operations table is first obtained so 
that the equilibrium of the ring as given by this solution vill provide 
a check on the whole setup. Joint G is considered fixed so that a 
unique solution to this system of equations is obtained} thus there are 
14 degrees of freedom to be considered. The 14 unknowns are found by 
the method of~ reference 9 to be: 


v A = - 605.73 X 10 “3 in. 

Wg = 40.825 x 10 “3 radian 
v 3 = 35 -lMt x 10 “3 in. 

Ug = - 300.06 x 10 "3 in. 

Wq = - 11.445 X 10 "3 radian 
v c = 664.55 x 10 "3 in. 

^ = - 72.282 x 10 “3 in. 

Wjj > - 22.975 x 10 "3 radian 
v D = - 94.734 x 10 '3 in. 

Up = 90.130 x 10 “ 3 in. 

Wg, = 6.2337 X 10 "3 radian 
v^ = - 42.621 x 10 ~3 i n . 

Ujj = 32.513 x 10 “3 i n . 
v F = - 44.648 x 10 “3 in. 


( 22 ) 


These displacements give the following values of the unknown moments 
and tangential reactions at A , F, and G cn the bars rather t hen on the 
Joints : 
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N A = - 3118.8 in. -lb 
T a = 402.51 lb 
Kg, = 105.43 in. -lb 
t f = - 182.57 11 
Nq. a - 371.06 in. -lb 
Rq. = 0.43 11 
T Q = 584.98 11 


> ( 23 ) 


Figure 9 is tie lending -moment diagram, for tie ring vitl these reactions 
applied. 

By examining the equilibrium of one -half the ring under these 
reactions and the external forces, the accuracy of the operations talle 
is established. Since R A and Bp are zero, the summation of forces 

in the vertical direction is simply: 


X'A = % = °*43 ^ 


( 24 ) 


The s umma tion of forces in the horizontal direction is: 


Z f h = t a - T F - T & = °* 10 115 


(25) 


The summation of moments about point G is: 
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X % » m a + t a ja 4 + — + 24(0. 70711 )j 

+ Wq. + % + Tj, ( 210(1 - 0.70711) - 1000(2 X 36 X 0.25882) (26) 

= - 3118.91 + 1*02.11.9(57.914.1) - 371.06 + 105. 1*3 
- 182.56(7.029!*.) - 18635.OI* 

= 17.45 in. -lb 


The equilibrium conditions for the half ring are approximately 
satisfied, the maximum percent error being a moment of less than 0.1 per- 
cent of the applied couple of 18,635 inch -pounds- It is considered that 
the accuracy of the operations table is established by this equilibrium 
check. 

Approximately 20 man-hours by an unskilled confuting -machine 
operator were required to solve this system of 14 equations. It is 
estimated that a skilled operator familiar with the Crout method would 
require about 10 man-hours . 

In applying the Crout method to this problem the coefficients of 
the linear equations are aBBumed to be mathematically exact and, there- 
fore, as many figures as could be carried on the 10 -bank computing 
machine are used throughout the computation. In thiB way an accurate 
solution Is obtained and the additional computing work is not great. 
Afterward the values of the unknowns can be rounded off to the physically 
correct number of significant figures. 

Use of the growing -unit method of solution cn. this ring is demon- 
strated as follows. This method Is described in detail on pages 39 to 46 
of reference 5 * It Is demonstrated on this new ring as an application of 
the procedure tcra ring with many intermediate Joints between the center 
line of symmetry points. In applying the growing -unit method to this 
ring the units are combined into bars of increasing length until dis- 
placements of all points are known such that the only unbalanced forces 
remaining act in the radial direction at A and I when unit radial dis- 
placements are undertaken at A and F. Then these forces at A and F can 
be eliminated by appropriate radial displacements of A and F and the 
final distorted shape determined. 

The first units to~be combined are AB and BC. In order to effect 
this combination, the displacements of B required to maintain the balance 
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of B during a unit radial displacement of A and unit radial, tangential, 
and rotational displacements of C must be determined. The displacements 
of B required to maintain the balance of B while point A is displaced 

radially 10 ”3 inch are given hy the equations: 

% = - 454.34*3 + 6.723873 - 78.ij.1iuB + 5.9020 X 10' 3 

% = 6.7238*3 - 12.093v B + 0.55690^ - 4.5778 x 10 " 3 

t b = - 78.411*3 + 0.5569OV3 - 84.510-Ug + 14.662 X 10 ~ 3 

The solution to these equations is: *3 = - 0.026434 X 10 “3 radian, 

V3 = - 0.38424 x 10 “ 3 inch, and Ug = 0.19549 X 10 “ 3 inch. 

If the forces and moments at points A and C due to a displacement 
v A = 10 ~3 inch and due to the foregoing displacements *3, 73, and u B 
are summed, the following equations are obtained: 



B A - - 2.66l8 lb 

Nq = IO.699 in. -lb 
Bq * - 1.8871 lb 
T c = 3.2614 lb 


(28) 


Ng, Rg, and T3 are zero since that is the condition satisfied by 
equations (27) . 

The displacements of B required to maintain balance at B during 
unit rotational, radial, or tangential displacements of C are determined 
in a similar manner and are collected in table 12. 

The forces and moments given in the last seven rows of this table 
constitute the influence coefficients for a new unit of the ring, namely, 
the segment ABC. This unit is not a bar, the center line of which Is an 
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arc of 1 a circle, "but rather one composed of two arcs of circles. This 
combining of units, extended until the entire ring is one segment, is 
the main principle of the growing -unit method. 

Each column, of table 12 represents a group displacement made up 
of individual displacements of points A, B, and C. Let these group 
displacements he identified by the Roman numeral given at the head of 
each column. For example, group II Is made up of the displacements 

wq = 10 " 3 radian, wg = - 0.21631 X 10 "3 radian, v-g = - O.2397I X 10 “3 inch, 

Ug = 0 . 7^197 X 10 “3 inch, and - Y c = u^> = 0. The moment at C, for 
instance, caused by the application of Xjj units of the group dis- 
placement II is then 


N c = - 389. 56X-J-,. 


(29) 


With a similar notation for all other forces and group displacements, 
equations (30) may be set up representing the requirements for equilibrium 
of joint C under the external forces acting at that -point, balance of B 
being maintained. 


N c = - 389.56XU 

R c = - 10.093Xj.j- 

T c = - 53 -77l3iii 


10.0932m - 53 - 771 *iv = 0 




8-75 2 9* iii 


10 • 615x T T t 


10.6153^. + 965-93 = 0 
5lf.l99x 1 y - 258.82 = 0 


► (30) 

J 


The solution to this system is Xjj = I.OL76, Xjjj = 217.61, and 
Xjy =* - 48.436, and the following forces and moments are introduced at 
A and D: 


R a = - 557-^3 lb 
= - 2666.6 in. -lb 
Rp = 184.66 lb 
Tp = 626.46 lb 




> (3D' 





MCA TN No. 1786 


39 


The forces and moments at D are added to the external forces applied to 
the ring at D and are balanced after the unit problem for the segment 
at ABCD is established. The R A force is not balanced until the complete 
ring is one segment and until the R A and Eg. residuals can be balanced 
together. 

The next unit to be considered is the combination of the ABC segment 

•with bar CD into the segment ABCD. The problem is to find the forces and 
moments at A, D, and E due (l) to a unit radial displacement of A with 
Joint D fixed and ( 2 ) to unit radial, tangential, and rotational dis- 
placements of D with A and E fixed. Joints B and C'are free to displace 
so as to maintain the balance at B and C in each of these four cases. 

By determining the magnitudes of xjj, xjjj, and xjy required to 
balance C in each of these four cases, the required displacements of 
both B and C are implicitly determined and the unit problem for segment 

ABCD solved* 


The magnitudes of the x TT , x T TT , and xjy operations required to 

balance Joint C when A is displaced radially 10 “8 inch and B permitted 
to displace so as to remain in balance are given by the following 
equations : 

Nq = - 389*58xjj - 10*0932:1x1 - 53 * 77 l z jv' + IO.699 = 0 ■"'j 


Rq = - 10 . 093x77 - 6 - 7529 z xii _ io.6i5xjy - 1.8871 = 0 > (32) 


T c = - 53 . 771 *n - 10 . 6 l 5 x irr - 54.199x3^ + 3-2615 =0 j 

The forces and moments at C to be balanced are given in group X in 
table 12 . The solution to these equations is Xjj = 0 . 021497 * 

Xjjj = - 0 . 53843 , and Xjy = 0 . 14430 . Use of these multiples of 

operations II, III, and IV and of a unit amount of group I results 
in the following forces and moments at A and D: 


R a = - 0.94505 lb 
— 6.7928 in* —lb 

► 

r d = - 0.74924.1b 
Td = 0.77749 lb 


( 33 ) 
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The forces and moments given by groups V, VI , VII, and VTII in 

table 13 . are the influence coefficients for segment ABCD. For example, 
the forces and moments introduced at A, D, and E due to a unit radial 
displacement of D with A and E fixed and with B and C in balance are 
given by VH. With these sets of coefficients it“is possible to 
balance Joint D while the balance of B and C Is preserved. The forces 
and moments to be balanced at D are (l) the external forces on the 
ring at D and (2) the forces and moments which are introduced at B by 
the balancing of C and which are given by equations (31) . The 
residuals to be balanced at~B are thus: 


Njj » - 2666.6 in. -lb 


% =■ - 965-93 + 18 ip. 66 = - 781.27 lb 


T D = - 258.82 + 626.46 = 367.64 lb 


“A 


( 3*0 


J 


The equations which condition the balancing of Joint B, from con- 
sideration of groups VI, VH, and VIII, are seen to be: 


Njj = - 346 • 88xy j - 16 . 697x^1 - 43 • 7 45 xy T TT - 2666 .6=0 


"A 


Rj-j — - 16 •697^'j “ 5 • - 12 * jj j- “ =0 ^ (35) 


T d = - 43. 7^52:71 - 12- 458xy JJ - 50 . 8l7x YIII + 367-64 = 0 


J 


The solution to these equations is Xyj = - 3.3IDO, Xyjj = - 328.09, 
and Xyjjj = 90.518, which give the following forces and moments: 


R a = 293.71 lb 
N e = 4889.4 in. -lb 
R e = - 522.69 lb 
Tg = - 149.69 lb 



( 36 ) 
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As in tlie "balancing of C, a tangential force and moment are introduced 
at A "by this "balancing of D, "but because of symmetry the equilibrium 
of A is not disturbed by these. The R^ forces will be balanced later 

and the residuals at E will be balanced when the influence coefficients 
for segment ABODE have been determined. 


In order to find the influence coefficients for bar ABODE, the 
forces and moments at A and E due to a radial displacement of A with E 
fixed and at A, E, and F due to unit radial, tangential, and rotational 
displacements of E with A and F fixed must be determined. By determining 
the magnitudes of groups VT, VII, and "VTII required to balance D in each 
of these four cases, the required displacements of B, C, and D and the 
required forces and moments are determined. 


The magnitudes of the groups VI, VII, and VIII required to balance D 

when joint A is moved radially 10 inch are given by the following 
equations : 


N D = - 346.88x^3- - l£. 697x^32 - 43-7 ^v ttt + 6.7928 = 0 

Rjj = - l6 . 6972^- -J- - 5* 6^00Xyj-p - 12 « ^t-^8xy-[--p-[- - 0-74924 = 0 


■"V 


> (37) 


T d = - 43.7453^-!- - 12. 4582^- - 9 Q» 8 lTVrTT + 0-77749 = 0 J 

The forces and moments at D to be balanced by groups VI, VII, and VIII 
are given by V in table 13 and are the constant terms in equation (37)- 
The solution of these equations is Xyj = 0.028042, Xyjj = - 0.42660, 
and Xy -j-p-p = 0.095744. The summation of forces and moments due to a 

unit magnitude of group V and the foregoing multiples of groups VT, VII, 
and VTII are: 


R a = - 0.36050 lb 
N e = 4.1055 in. -lb 
% = - 0.31703 It 
T e = 0.19226 lb 




> (38) 
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In a similar manner the complete set of influence coefficients for 

segment ABODE is determined and is given in table 14 . For example, the 
forces and moments in group 21 are the forces and moments introduced 
at— A, E, F, and G by a unit radial displacement of E with A, F, and G 
fixed and with points B, C, and D free to displace so as to~ remain in 
equilibrium. With these influence coefficients Joint E can be balanced 
while the balance of B, C, and D 1 b preserved. 

The forces and moments to be balanced at Joint E are those introduced 
by the balancing of Joint D with groups VI, VII, and Till and are given 
by equation (36) . 


The equations in x^., x^, and x-jq-j balancing Joint E under these 
loads are: 


= - 533-92* z - 38-0 99x53- 

Rg = - 38.099 x x - 49.295x33- 
T e = - 29.579x3- - 53.432x33- 


29^79x XII + 4889.4 
53.432x xri - 522.69 
76.322 x- rrT - 149.69 



0 > ( 39 ) 

0 


The solution to these equations is x^- * 11 . 184 , x-jg- = - 51.518, and 
Xjji = 29.771 and the forces at A, F, and G introduced by this 
balancing of E are: 


% = 67.974 lb 
Rp = - 266.68 lb 
% = 70.668 lb 



( 40 ) 


The tangential forces and the moments introduced at A, F, and G are not 
considered in this balancing of the half ring, since these are 
equilibrated by the forces and moments from the other half of the ring. 

The final combination of units will be the combination of bar EF 

with the unit ABODE. When this union is effected, the influence coef- 
ficients for the half ring as a unit will have been determined and the 
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radial forces at A and F can "be balanced s imultaneously . The radial 
forces at joints A and F due to a unit radial displacement of A with F 
fixed and to a unit radial displacement of F with A fixed must he 
determined. In both cases joints B, C, D, and E are displaced so as to 
remain balanced. 

The equations giving the magnitudes' of groups X, XI, and XII 
required to balance Joints B , C , D, and E when Joint A is displaced 
radially as in group IX are: 


% = - 533 - 92 x t - 38.0992:^ - 29. 579x yTT + 4.1055 = 0 


txe 


XII 


% = - 38.099x z - 49.295x2! - 53.432x211 - 0.31703 = 0 V (4i) 


L E 


= - 29.579x2 - 53*^32x21 ~ 76.322x2!! + 0.19226 = 0 


The solution to these equations is X2 = O.OO94398, Xyj = - 0 . 05 l 8 o 4 , and 
X2!! = 0.035128 and the forces introduced by a unit magnitude of IX and 
by these multiples of groups X, XI, and XII are: 


R A = - O.29857 lb 
RF = - 0.33361 lb 
% = 0.035176 11 


► ( 42 ) 

.J 


The equations giving the magnitudes of groups X, XI, and XII required 
to balance Joints B, C, D, and E when Joint F is displaced radially 

1D~3 inch are: 


Ng — - 533*92X2- - 38.O99X2! ” 29-579X2II 


— - 38.099x2 _ 49 * 295x2! - 53 • 432 xy-,- t 
T e = " 2 9 * 579 x z - 53 *^ 32 x 22 - 76.322x211 


5.9020 = 0 
4.5778 = 0 


-\ 

► ( 43 ) 


14.662 = 0 
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The solution to these equations is xj = - O.OI 7182 , = 0-50354, and 

x^-p = - 0-53797 and- "the forces introduced atr^A, F, and G by a radial 

displacement of F of 10 ”3 inch and by the foregoing multiples of 

groups X, XI, and XII are: 


R A = - 0.33361 Tb 
% = - 1.44-70 lb 
Rq. = 1.1156 lb 


"A 

> (44) 


The forces given by equations (42) and (44) represent the influence 
coefficients for the entire half ring and are labeled groups XIII and 
XlV, respectively. These forces permit calculation of the multiples of 
groups XIII and XIV required to balance the radial forces at A and F- 
These forces are the total forces remaining from the balancing of C. D, 
and Ej R A is given by the sum of the R A forces of equations (31), 

( 36 ), and (40) and is: 

R A = - 557*43 + 293.71 + 67.974 = - 195*75 lb 


The Rg. force is the force introduced by the balancing of E alone 
and Is given by equation (40). It is: 


R F = - 266.68 lb 


The equations giving the magnitudes of groups XIII and XIV required 
to balance joints A and F under those loads are: 


S A 

Rj. 


0.29857x22!! 


0 . 3336l X7-|"rT 


0 . 33361 x 2 ^ - 195*75 = 0 
l.4470xjxv “ 266.68 = 0 i 


( 45 ) 


The solution to these equations is ^xill = ” 605*73 and 

44.646. The radial force at G introduced by this balancing 


x XXV 
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is - 71-114 pounds, Nut the Bq, force given by equation ( 4 o) in the 

"balancing of E is 70.668 pounds . The difference "between the two, 

- 0*446 pound, is considered negligibly small compared to the applied 
loads of 1000 pounds • 

With the balancing of joints A and E and the substantiation of the 
balance at G-, the entire half ring, is balanced. The total deflectiQns 
in each degree of freedom can now be calculated and used to determine 
the unknown bending moments and tangential forces at A, F, and G. In 
order to calculate tfiese deflections the balancing equations (30), (35), 
(39), and ( 45 ) give the magnitudes of the group operations involved 
while the equations determining the group influence coefficients give 
the individual operations involved in each group. 

Table 15 gives the magnitude of all group displacements from I 
to XIV implied in a unit application of any one group. For example, 
row X in this table indicates that a unit magnitude of group X (that is, 
Xj- = 1) is equivalent to the sum of the effects of Xyjjq = 2. 1885, 

Xy TT = - 4.6760, Xy.j. = - 0.l6l86, and w^ = 10 “3 radian, or the sum of 

the effects of % = 2-9219, % = 1.7552, x XI = - 0.19736, 

w E = 10 "3 radian, w-p = - O.16186 x 10 "3 radian, Vjj = - 4.6760 x 10 3 inch, 

and u D = 2.1885 x 10 “ 3 inch- During the solution of the problem the 
magnitude of group X which was explicitly used was 11 . 184 , as given in 
the last column of table 15* 

From table 15 the total magnitudes of each group operation may be 
found. For example, the total magnitude of group VI is: 

x V i = (1)(- 3-3100) + ( 0 . 028042 ) (0) + (- 0.16186) ( 11 . 184 ) 

+ (0.037474) (- 51.518) + (- 0.0022427) (29.771) 

+ (0.024494) (- 605.73) + (0.022857) (- 44 . 646 ) ( 46 ) 

= - 22.975 


The total displacement wp is: 


w D = ( 2 yi) X 10 “3 = - 22.975 x 10 "3 radian ( 46 a) 

Similarly the displacements of all points except point B may be 
calculated from table 15 and are given in the last row of that table- 
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Point B was displaced during the application of groups I, II, III, 
and IV, and therefore the magnitude of its displacement must he calcu- 
lated as indicated in the following example: 

Wg = (- 0.026434 -x 10"3)(- 605.73) + (- 0.21631 x 10“3)(- 11.444) 

* + (0.035554 x 10 ”3) (664.55) + (0.017847 x 10 "3) ( _ 72.286) (47) 

= 40.825 X 10 “3 inch 


where the first" number in each product is the magnitude of wg involved 
in each unit application of groups I, II, HI,* and IT, respectively. 

The total ' displacements used are assembled in equations (47a). 

v A = - 605.73 x 10 "3 in. 


Vg = 40.825 X 
v B = 35-144 x 

Tig = - 300.06 

Vp = - 11.444 
Vq = 664*55 X 

Uq = - 72.286 

w D = - 22.975 
v D = - 9^-731 
Ujj = 90.127 X 
Wg = 6.2331 X 
Vg = - 42.620 

Tig = 32 . 5 H X 

yg = - 44.646 


10 “3 radian 
10 "3 in. 
x 10"3 in. 

X 10 “3 radian 
10 “3 in- 
X 10 "3 in. 
x 10 “3 radian 

X 10 “3 i n . 

IO "3 in. 

10 “3 radian 
X 10 "3 in. 

10 "3 in. 
x 10 "3 .in. 


> (47a) 
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These total displacements constitute the unknowns of the system of 
equations given hy the operations table and the external forcesj 
comparison between this growing-unit and the matrix calculus solutions 
given by equations (Vfa) and (22), respectively, indicates good agreement 
for the displacements. In fact, the forces and moments given by the two 
methods differ by less than 1 percent and therefore are given only for 
the matrix method (equation (23)). 

Several general remarks are made about the grow ing -unit method: 

(a) In determining the influence coefficients and in balancing the 

external forces and moments, sets of equations with the same 
left-hand sides but with different constant terms are used 
several times. This simplifies solution of the equations 
and reduces the computational work considerably. 

(b) In order to obtain sufficient accuracy of solution for rings 

with many Joints, calculating machines must be usedj five 
significant figures were carried throughout the calcu- 
lations. However, on the simpler rings such as the circular 
ring and the egg-shaped ring discussed previously, slide - 
rule accuracy for determining the displacements in a 
combined operation is probably sufficient for engineering 
purposes • 

(c) A check on the influence coefficients for composite bars is 

obtained by applying Maxwell 's theorem of reciprocal 
deflections. This Is a valuable device for assuring 
accuracy at each stage- 

In applying the new relaxation procedures to this ring, it would 
have been possible to. use the general method described for the egg- 
shaped ring, that Is, to consider the residuals after each operation and 
develop a satisfactory combined operation to reduce as many residuals as 
possible. However, the number of degrees of freedom involved in this 
ring is large and, therefore, the number of residuals to be considered 
in testing the efficacy of a particular operation is large. 

The loading on the ring provides a clue to overcoming this difficulty. 
No external loads are applied at A, B, E, F, and G$ moreover. A, F, and G 
are points along the center line of symmetry. Therefore, if in balancing 
D and E the balance at the other Joints is preserved by suitable 
displacements, attention is fixed on the two Joints D and E and the pro- 
cedure described for the egg-shaped ring can be used effectively. It 
will be recognized that this procedure is essentially a combination of 
growing-unit and relaxation methods of solution. 

In executing the proposed method the bar ABCD, free only to displace 
radially at A and fixed at D, is considered first. The equations giving 
the displacements of A and B required to maintain balance of these points 

while Joint C Is rotated through 10"3 radian are: 



B a = - T-1310t a + 5-9020^ 


4.577&tb + 14-6621^ = o 




1 % = 5 * 9020 t a - lt.54.3lj.Wg + 6.723873 - 78.411113 - 38.489 X 10 "3 = 0 

% = - ^- 57 TSv a + 6.7238^ - 12.093t b + 0.55690^ - 1.8576 x 10~ 3 = 0 

Eg = 14 . 662 v a - 78.4llwg + 0.55690^ - 84 . 5 l 0 ug + 45.876 x 10 _3 . 0 


> m 




The solution to these equations is v A - 4.0195 X 10 “ 3 inch, wg = - 0-32255 x 10 ' 3 radian, 

V B “ “ 1‘7 3 ^2 x 10 3 inch, and Up = 1.5277 X 10 “ 3 inch. The Be displacements combined with the 
unit rotation of C yield: 


N c * - 3^- 56 in. -lb 

Bq = - 17.678 lb 

T c = - 1(0-663 1T> 

Up = - 38-489 In. -lb 

Rp = 1.8576 lb 
Tp » 45-876 lb 


> (49) 


J 


The moment and tangential force introduced at A are not considered until the balancing of the 
ring la complete. 


£ 


l 
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In a similar manner the forces and moments for unit radial and 
tangential displacements of C are determined, as shown in table 1 6. The 
forces and moments given by groups XV, XVI, and XVII constitute the 
influence coefficients for the displacements of C with A and D fixed and 
with joints A and B balanced. Use of these coefficients permits focusing 
of attention on C and D, the joints at which the external forces are 
applied when C is being balanced. 

The forces and moments introduced at C and D when D is displaced a 
unit amount in each degree of freedom and when E, E, and G are displaced 
so as to maintain the balance thereof are calculated and shown in 
table 17 • 

Table 18 is an operations table consisting of unit magnitudes of 
group operations XV to ZX. Table 19 is -the relaxation table for this 
ring which uses these group operations. The external forces applied 
at C and D are given in the first row of table 19 - 

A discussion of each step in the relaxation process Is given as 
follows . 

Step 1 . - Because of the antisymmetry of the loading and of the 
quasisymmetry of the ring about a horizontal axis operations = 1 

and x jjjr = -1 are applied as a first approximation to the deflected 

shape* The forces and moments introduced are as given In the following 
table : 


Eorces and 

" — - mrnnaritH 
Operat ion""-^^^ 

N 0 


T c 

n d 

=0 

*D 

SJ 

H 

II 

I- 1 

-17.678 

-5.4150 

-12.928 

-1.8576 

-2.2277 

I3.78I 

(XIX) = -1 

-I.8576 

2.2277 

13.781 

- 11-535 

6.2441 

-12.706 

Operation K = 1 

-19.536 

-3-1873 

0.853 

- 13-393 

4.0164 

I.O75 


Operation K is used to balance the residualj the same operation 

reduces the other force residuals but introduces large Nq and 
Np residuals. 

Step 2 . - In order to reduce these moment residuals an ant isymmetr leal 
combination of Wq and w D is made, as shown as group operation L: 
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Forces and 
'''--^mDments 
Operat Icri~' — 

% 

BC 

T C 

% 

e d 


(W) = 1 

- 346.56 

-17.678 

- 4 o . 662 

-38.489 

1.8576 

45.876 

(XVIII ) = 1 

-38 . 489 

-I.8576 

45.876 

-392.46 

H -535 

- 42.305 

--> Operation I> = 1 

-385.04 

-19.536 

5.214 

-430.95 

13-393 

3-571 


However, use of operation L by itself would reintroduce large Rq and 
Rjj residuals, and therefore operations E and L are combined so that - the 
Eg residual will he smaller and the residual eli m inated, as shown 

in the following tables 


Forces and 
'--^moments 

Operation 

*0 


T c 


=0 

T 

d 

Operation L = 1 

-385.04 

- 19-536 

5-214 

-430.95 

13.393 

3-571 

- 3-3346 x Operation K 

65.145 

ID. 628 

- 2.844 

44.660 

-13.393 

-3.5846 

— > Operation M = 1 

- 319-90 

-8.908 

2.370 

-386.29 

0 

- 0.0137 


The new force residuals introduced by operation M are less than 30 percent 
of the original residuals and, therefore, the rate of convergence is felt 
to be adequate. 

Step 3 * - The radial residuals at C and D have the same sign and, 
therefore, symmetrical displacements vq and v D are undertaken- It is 

seen that such a combination would introduce large tangential residuals 
at C and D. Therefore, a tangential displacement of C (D could have been 
chosen Instead) such as to eliminate the Tq and Tj> forces is under- 
taken. The forces and moments introduced by the individual displacements 
and by the combination are denoted as operation N. 
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Forces and 
^"^moments ' 
Operation^\^^ 

N c 

R c 

T 

C 

n d 

*D 

T 

D 

(XVI) = 1 

-17.678 

-5.i4.150 

-12.928 

-1.8576 


13.781 

H, 

II 

H 

I.8576 

-2.2277 

- 13-781 

11-535 

-6.2kki 

12.706 

-0.53201 x (xvii) 

21.632 

6.8773 

26.709 

-2k. J W6 

7-3316 

-26.587 

• ■> Operation IT = 1 

5.8U 

-0.7653 

0 

-ik. 730 

-l.lkoo 

-0.100 


The use of operation N reduces substantially all the residuals except 

Step 1 - - In order to reduce Kq and at the same time keep the T c 
and Tj) residuals small, a combination of groups XV and XX is made. 
Group XX is included since a force increasing the residual Tj> vould he 
introduced hy the use ‘of XV alone • 


Forces and 

— TnnTnA-ntfl 
Operation" — 

N C 

*C 

T c 

% 

*1 

T D 

(zv) = 1 

-3k6 . 55 

-17.678 

-ko.662 

-38 . k89 

1.8576 

k 5-876 

0.91279 X (XI) 

kl.875 

12.579 

k5 . 616 

-38.616 

11.597 

-k5 • 876 

Operation 0=1 

-30k. 67 

- 5.099 

k- 95 k 

-77. 1 D 5 

13 -k 55 

0 


Steps 5 and 6 . - After operation 0 is used, the largest force residual 
is fippm y 1 nw.te 1 y 6 percent of the applied forces and the moment residuals 
are small. It was considered desirable to reduce further the force 
residuals. Therefore, operation I was used again so as to reduce E^ 

the largest force residual, and then XVII was used so as to reduce the 
result ing T c residual. After this sixth step the largest residual of 

k percent of the external force is considered small enough. 

A check table using the total displacements is used as a check on 
the accuracy of the combined operations and on the relaxation table. The 
total individual displacements are calculated as discussed in the previous 
two examples and are as follows: 
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v A = - 596.18 x 10 " 3 in. 
v B = 36.779 X 10 “3 radian 
7 b = O.22394 x 10 "3 in. 

■Ug = - 304.69 x 10 "3 in. 
v q - ~ 1^*32 X 10 “3 radian 
Yq =s 561.66 x 10 ~3 
Uq = - 114.61 x 10 "3 in. 
w D = - 18.4 X 10 ~3 Indian 
v D = - 133.66 x 10 “3 in. 
u D = 3.7242 x 10 “3 in. 
v E = 7.4813 x 10 "3 radian 
v E = 3O.5O7 X 10 "3 in. 

Ug = - 41.708 x 10 '3 in. 

= 59*979 x 10 "3 in. 
v G = 114.52 X 10 “3 in. 




> (50) 




It is pointed out that certain of these displacements differ 
considerably from those given by the exact solutions of the matrix 
calculus and growing -unit methods, mainly “because the re lax ation solution 
is approximate and in it Joint G- is permitted to^displace radially. 

The unknown reactions given by the foregoing relaxation procedure 

are: 
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W A = - 2851.1 In. -lb 
t a = 380.21 11 

Up = lJ»0.l6 In. -11 
Tj, = - 224.09 lb 
W Q = - 440.45 in. -11 
Tq. = 648.18 11 


> ( 51 ) 




Consideration of the equilibrium of the half ring gives: 


F g = 380.21 + 224.09 - 648.18 = - 43.88 11 

51 % - 0 



2851.1 + 380.21(57.941) - 440.45 + HiO.16 

» 


> ( 52 ) 


- 224.09(7.0294) - 18,635 


= -1331.8 in. -11 


The moment equilibrium unbalance is approximately 7 percent of the 
applied moment and is considered satisfactory for engineering purposes. 
If a more accurate representation of the final deflected shape and 
consequently of the lending -moment diagram is desired, several more 
operations in the relaxation talle could le undertaken and the residuals 
at C and D further reduced. 

The lending -moment diagram given ly the reactions of equation (51) 
is shown in figure 9 along with that of the exact solutions. The 
external unbalanced moment of 1331.8 inch-pounds is applied linearly 
along the ring as a distributed moment. If this unbalance is not 
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distributed in this manner, it would "be concentrated at- either Joint A 
or Joint F, depending on the direction in which the bending moments are 
calculated, and would lead to large errors in the bending moment in the 
neighborhood of that Joint. It is seen from figure 9 that the agreement 
between the exact and relaxation solutions is good. 

It is pointed out that^— by slightly modifying the determination of 
the influence coefficient for Joint ! when E is fixed and F and G free 
to displace radially, a table similar to table 17 could be established 
and solved by matrix calculus methods. The slight modification is to 
walfft V G = 0 in the equations corresponding to table 17. Such a 

solution is essentially the growing-unit method, except that the ring Is 
combined from Joints C and D to A and F, respectively, rather than from 
A to IF • The total displacements in each degree of freedom, will be the 
same in each approach. 


CONCESSIONS 


This report contains recommendations as to the choice of the most- 
expeditious method of solution of the s Imultaneoub linear equations 
represented by the operations table and the external loads. The 
operations table is first established in accordance with Southwell 's 
suggestions and, together with the external loads, defines completely 
the problem of stress distribution in a reinforced panel or of the 
moment distribution in a fuselage ring. However, the following 
generalized suggestions can be made: 

1. In most reinforced panel problems the use of the relaxation 
procedure Is advantageous. 

2 . Solution of the equations defining a reinforced panel problem 
by thabtis of the electric analogue is advisable when many closely related 
problems have to be investigated. 

3. Ring problems are best solved by matrix methods. 

4 . In very complicated ring problems a combination of matrix methods 
with the growing -unit and relaxation methods may become advisable. 


Polytechnic Institute of Brooklyn 
Brooklyn, N. Y., June 25, 1947 
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TABLE 2.- RELAXATION TABLE FOE REINFORCED PANEL - PROCEDURE 1 

[Cycles of operations shown should "be repeated until residuals are 
considered negligibly small. Forces are in Ibj 

displacements, in in. x 10“^] 


— ^Force 
Op eratioSr^--^^ 

Q 

*B 

y e 

H 

□ 

D 

% 

*0 

Erteraal forces 
T Bloc 3 c 1 " " 2 ‘ 5 

-120 

10 

-10 

20 

-20 

20 

-20 

60 

10 

. 60 
-10 


-no 

-10 

20 

-20 

20 

-20 

70 

50 

t e = 2.35 

no 

5 

-238 

9 

no 

5 

0 

0 


0 

-5 

-218 

-n 

130 

-15 

70 

50 

Tj = 4.65 

0 

0 

218 

9 

-473 

19 

218 

9 


0 

-5 

0 

-2 

1 

■n 

288 

59 

* 7-33 

0 

0 

0 

0 


mm 

-372 

14 


0 

mm 

0 

-2 

0 

18 

-84 

73 

T BlocE 2 = 3 *5 

14 

El 

28 

-28 

28 

-28 

14 

-14 


■9 

-19 

28 

-30 

28 

-10 

-70 

59 

▼j. = 0-371 

M 

19 

1 

-ki 

1 

19 

0 

0 

• 

15 

0 

29 

-71 

29 

9 

-70 

59 

T K = 1.385 

0 

0 

3 

71 

5 

-153 

3 

71 


15 

0 

32 

0 

34 

-144 

-67 

130 

Tq = 2.81 

0 

0 

0 

0 

5 

144 

5 

-156 


15 

0 

32 

0 

39 

0 

-62 

-26 

T Block 1 = 1,0 

-4 

k 

-8 

8 

-8 

8 

-4 

4 


n 

k 

2k 

8 

31 

8 

-66 

-22 

t e = - 0.235 

-n 

0 

2k 

-1 

-n 

0 

0 

0 


0 

k 

k8 

7 

20 

8 

-66 

-22 

Tj = -1.025 

0 

0 

-k8 

-2 

io4 

-4 

-48 

-2 


0 

n 

0 

5 

124 

4 

-114 

-24 

v H . -2.65 

0 

■a 

0 

0 

-124 

-5 

135 

-5 


0 


0 

5 

0 

m 

21 

-29 
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TABLE 3 . - RELAXATION TABLE FOR REINFORCED PANEL - PROCEDURE 2 




[Forces are In lb} displacements. In In. x 10 


• 4 ] 


^"‘"‘—'-^.^Force 

Operatiao^-^^^ 

Y A 

y b 

y e 

y f 

Y J 

y k 

y n 

Y o 

External forces 
v Block 1 = -2 ' 5 

-120 

10 

-10 

20 

-20 

20 

-20 

60 

10 

60 

-10 

v A =» -2.165 

-no 

110 

-10 

-4.3 

20 

-101.3 

-20 

-4.3 

20 

0 

-20 

0 

70 

0 

50 

0 

(1) = -1.482 

0 

5*9 

-14.3 

- 5*9 

-81.3 

81.3 

-24.3 

-8.9 

20 

-69.4 

-20 

- 3*0 

70 

0 

. . J 

50 

0 

(2) = -0.901 

5*9 

3*6 

-20.2 

- 3*6 

0 

7.2 

- 33*2 

- 7-2 

- 49.4 

49.4 

-23.O 

- 5-4 

70 

- 42.2 

50 

-1.8 

v Block 2 = " 1,8 5 ° 

9*5 

-7.4 

-23.8 

7-4 

7.2 

- 14.8 

- 40.4 

14.8 

0 

- 14.8 

- 28.4 

14.8 

27.8 

- 7-4 

48.2 

7.4 

7b = -o - 2 97 

2.1 
“0 • 6 

- 16.4 

16.4 

-7.6 

-0.6 

-25.6 

-15.2 

- 14.8 

0 

- 13*6 

0 

20 . 4 - 

0 

55 *6 
0 

(3) ■ -0.689 

1*5 

-2.8 

0 

2.8 

-8.2 

- 4.1 

- 40.8 

40.8 

- 14.8 

-1.4 

- 13*6 

- 35*3 

20.4 

0 

55*6 

0 

( 4 ) = -0.825 

- 1*3 

- 3*3 

2.8 

3-3 

-12.3 
—6 . 6 

0 

6.6 

-16.2 

-5.0 

-48.9 

48.9 

m 

55.6 

-42.3 

tr„, , , a -I.08 
Block 1 

- 4.6 

4.3 

6.1 

- 4.3 

-18.9 

8.6 

6.6 

-8.6 

-21.2 

8.6 

0 

-8.6 

BE 

IBl 

ESI 

(1) = -0.188 

-0.3 

0.8 

1.8 

-0.8 

-10.3 

10-3 

-2.0 

-1.1 

-12.6 

-8.8 

-8.6 

- 0.4 

23.1 

0 

9.0 

0 

( 2 ) = - 0.390 

0-5 

1.6 

1.0 

-1.6 

0 

3*1 

- 3-1 

- 3-1 

- 21.4 

21.4 

- 9-0 

-2.3 

23.1 

-18.3 

9.0 

-0.8 

T Block 2 = 

2.1 

-1.6 

-0.6 

1.6 

3*1 

- 3*3 

-6.2 

3*3 

0 

- 3*3 

-n.3 

3*3 

ayfei 

m 

8.2 

1.6 

( 4 ) = -0.135 

0.5 

-0.5 

1.0 

0.5 

mm 

Baa 

-2.9 

1.1 

- 3*3 

-0.8 

- 8.0 

8.0 

3.2 

-0.3 

6.6 

-6.9 

> 


IB 

KB 

-1.8 


0 

2.9 

-0.3 


V --*. V' 
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TABUS 4 - - RELAXATIOB TABUS FOR REINFORCED PANEL - FIXED ENDS 
[^Forces are in Ibj displacements, .in in. x 10 ~^j 


^~^''~'-'~^Force 

Operaticn^~-~^_^ 

External forces 

y a = -2.36 

(1) 

= -2.02 

(2) 

= - 1.723 

V B 

= - 0-357 

( 3 ) 

= -0.83 

( 4 ) 

= - 0-959 

V A 

= 0.1375 

(l) 

= 0.124 

(2) 

= -0.0292 

V B : 

= 0.143 

( 3 ) 

= 0.265 

( 4 ) 

= 0.228 



-120 

120 - 4*7 




8.1 - 12.8 
6-9 - 6-9 


15.0 -I9.7 
-0.8 19.7 



0 

13-8 


13-8 

- 0.8 


-16.8 

-13.8 


.94. 4 - 4.0 

94 . 4 -10.3 



-80.6 


-80.6 


14.3 -80.6 

42.3 




2.1 

1.8 

0 

-1.8 

3-9 

-1.8 



0 

-82.5 

-52.5 

0 

0.2 

-82.5 

-52.5 

0.2 

-0.3 

-82.5 

-1.4 

-52.5 

-0.1 

-0.1 

-83.9 

-52.6 

-0.1 

13.6 

-83.9 

-52.6 

13-5 

- 13-5 

-83.9 

0.5 

-52.6 

11.7 

0 

-83.4 

-40.9 


^NACA 


















































































TABLE 5.- OPERATIONS TABLE FOR REINFORCED PANEL - GROWING-HHIT METHOD 
[Forces are In Ibj displacements, in in. x 10 - 4j 



Forces and 

Operation 

— — ^jKJ^nbs 





(9) 0.0362 x (2) 

(10) (1) + (9) 

(11) O.921 x (1) 

(12) 0.0695 X (2) 

(13) (3) + (11) + (12) 
(llf) 0.0758 x (1) 

(15) 0.923 X (2) 

(16) ,(4) + (llf) + (15) 

(17) 0.828 x (13) 

(18) 0.158 x (16) 

(19) (5) + (17) + (18) 

(20) 0.170 X (13) 

(21) 0-845 x (16) 

( 22 ) ( 6 ) + ( 20 ) + ( 21 ) 

(23) O.79 x (19) 

(24) 0.19 X (22) 

(25) (7) t (23) + (2lf) 

(26) 0.21'X (19) 

(27) 0.81 X (22) 

(28) (8) + (26) + (27) 


0.07 

- 50.7 

-U6.9 

0.1 

0 

-3-9 

I.85 

0 



46.8 

2 

- 101.6 

4 

46.8 

2 


O.O7 

46-9 

4-3.2 

0.1 

-58.3 

3-5 

1.85 

9-4 

-48.3 

1-5 
0 ; 
-9.9 

7-9 

0 


51.2 

4 

-110.4 

2 

51.2 


1.86 

3-8 

1.85 

3-5 

9-4 

0.2 

47.6 

-62.6 

7-9 

-9.9 

0 

1.6 

-52.8 

0 


46.8 

2 

-101.6 

4 


46-8 

2 

2 

51.2 

37-0 

1-59 

0.38 

9.83 

-64-2 

15.4 

9-8 

0.ll2 

1.6l 

41.2 

15-4 

-68.8 
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TABLE 6. - RELAXATION TABLE TOE KEIKFOECED PANEL - GSBOWUTG-ENIT METHOD 


m * 


J 


Forces are In lbj displacements, in In- x 



i 


Force 

Operat lon\^_ 

Y A 


y e 

b 


% 

% 

Y 0 

y e 

Y s 

External loads 

-2.37 X (10) 

-120 

120 


-111 

-9 







-1.975 x ( 13 ) 
-o.W* X (16) 

0 


-111 

115.2 

4.2 

-9 

-18.8 

-92.5 

-0.9 

4-0 

-22.7 





-1.657 X (19) 
-0.71*2 X (28) 



0 

0 

-934 

123.7 

-10.2 

-26.7 

-22-9 

1*9-6 

-77.6 

-1.5 

-3*3 

-37-9 



-1,1*5 x (25) 
-0.925 x (28) 





0 

0 

-79.1 

93.2 

-ll*.2 

41.2 

-224 

63*6 

47.9 

-1.8 

-2.9 

474 








0 

0 

„ 

-69.7 

-50.3 


o\ 

H 
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TABLE 7. - OPERATIONS TABLE FOR CIRCOLAR RUKJ 


Forces and 
^\^Bsoiaents 

Operat Ion'"'-\__ 

*A 

(in. -11) 

t a 

(11) 

% 

(in. -11) 

% 

(11) 

% 

(11) 

*0 

(in. -11) 

T C 

( 31 ) 

(I) v A = ID "3 radian 

-281.95 

49.079 

-29.966 

4.733 

64.675 


■ 

(2) u A = ID "3 in. 

-49.079 

-52*296 

64.675 

- 22.441 

51*516 



(3) Wg «= 10 “3 radian 

-29.966 

64.675 

439*849 

31*443 

- 50.642 

56.5H7 

6.632 

( 4 ) * 10 “3 in. 

- 4-733 

- 22.441 

31 * 4-3 

-12.338 

20.14 

8.842 

0.524 

(5) Ug = 10 "3 in. 

64.675 

51.516 

- 50*642 

20.14 

-52.618 

6-632 

O.0685 

(6) v c = 10 "3 radian 



56 * 511 ? 

8.842 

6.632 

-157.899 

-I.563 

(7) = 30 "3 in. 



6.632 

0.524 

0.0685 

-I.563 

-0.322 




4 


J 


f 


TABUS 8. - RELAXATION TABLE FOR CIRCULAR KING- 


Forces and - 
^\moinent8 

Operation's. 

*A 

(In- -lb) 

t a 

(lb) 

% 

(In. -Ib) 

% 

(») 

t b 

(lb) 

— 

*0 

(in. -lb) 

T c 

(lb) 

External Forces 

- 1.84 

-8.75 

- 55*0 

59-5 

38-1 

- 53-1 

-23.9 

-O.OO778 x (A) 

1.84 

0 

O.7 

-0.1 

-0.1 

0 

0 


0 

- 8*75 

- 5* -3 

59-4 

38.0 

-53.I 

- 23-9 

-0.2 x (B) 

0 

8.75 

- 14.0 

4-3 

-8.0 

0 

0 


0 

0 

-68.3 

63-7 

30.0 

- 53-1 

- 23-9 

- 0-353 x (C) 

0 

1 0 

-8.6 

-2.2 

-2.2 

53-1 

0 


0 

0 

-76.9 

6I.5 

27.8 

0 

- 23*9 

-77.8 x (D) 

0 

0 

-472 

- 34.0 

-0.2 

0 

23-9 

-2.98 x (E) 

0 

0 

-549 

27.5 

27.6 

0 

0 

0 

0 

549 

-27.5 

-27.5 

0 

0 


0 

0 

0 

0 

0.1 

0 

0 

Check-table 

results 

• 

0.0171 

0.0030 

O.5095 

- 0.0043 

- 0.0053 

-0.0726 



-0.0121 


a\ 

LO 
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TABLE 9. - GEKR£tIQff8 TABLE TOE &3G-SHAPED BUR? 


N. Forces and 
>s v ucjnents 

Operation 

B A 

(lb) 

% 

(in. -lb) 

% 

(lb) 

% 

(ID) 

«C 

(In. -lb) 

*0 

(lb) 

T c 

(lb) 

% 

(lb) 

( 1 ) t a » in^ In. 

-3.34833 

8.92216 

-2.6961* 

3-96771 





(2) a 1D"* ; rail an 

8.92216 

-327.866 

11-4697 

- 13.1014 

- 6 l . 242 

8.10267 

0 


( 3 ) Tb = In. 

-2.6o6l| 

11-4697 

-4.00991 

3-4352 

-8.10267 

O.66158 

0 


00 ng * lfl' 1 | In. 

3-98771 

- 13-2014 

3-4352 

- 30.9566 

0 

0 

26.2058 


(5) Ho m 10 ”^ radian 


- 61.242 

-8.10267 

0 

-288.367 

-2.95622 

-5.24667 

- 7 - 352 ^ 

(6) r c « lo 4 m. 


8.10267 

0.66158 

0 

-2.95622 

-1.90205 

-0.88929 

- 1.11900 

(7) ^ - 10 in. 


0 

0 

26.2058 

- 5.24667 

-0.88929 

-27-0833 

- 1-0375 

(8) v D - ao 4 in. 





-T -3524 

-1.31900 

-1.0375 

- 1.2400 
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TABLE 3 D.- RELAXATION TABLE FOR EGG-SHAPED RING 


N. Forces and 
^\moinents 

Operation^. 

• 

b a 

W 

% 

(In. -lb) 

% 

(31) 

% 

(lb) 

. 

% 

(In. -lb) 

*0 

(lb) 

T c 

(lb) 

r d 

(lb) 

External forces 

-500 

0 

0 

0 

0 

0 

0 

-500 

- 149.2 x (1) 

500 

-1330 

402 

-592 

0 

0 

0 

0 


0 

-1330 

402 

-592 

0 

0 

0 

-500 

-403 x (8) 

0 

0 

0 

0 

2960 

451 

4 l 8 

500 


0 

-1330 

402 

-592 

2960 

451 

418 

0 

1336 X (F) 

0 

0 

.0 

29 

-2485 

-451 

63 

0 


0 

-1330 

402 

-563 

475 

0 

48 i 

0 

256 x (I) 

0 

1615 

-402 

62 

0 

138 

-8 

0 


0 

285 

0 

-501 

475 

138 

473 

0 

409 X (F) 

0 

0 

0 

9 

-760 

-138 

19 

0 


0 

285 

0 

492 

-285 

0 

492 

0 

-9.4 X (J) 

0 

24 

-2 

493 

9 

0 

-492 

0 


0 

309 

-2 

1 

-276 

0 

0 

0 

Check table 

-0.465 

309.971 



-2.231 

0.938 

-280.405 

-0.202 

- 0.351 

-0.279 




TVCHTI 11. - OmtXZLOB TAB 13 1UR OfJOrSEABSD EUK 


Joroee and 



(») 

(to. - 31 ) 

(11) 

( 11 ) 

(to. - 11 ) 

(U>) 

( 11 ) 

(to. - 11 ) 

Tj^ *S 10 ^3 12 X # 

-7.13m 

5-9020 

4 - 57(8 

31 . 66 s 





njj - U 0"3 rallan 

5.9000 

-W- 3 + 

6-7238 

-78.IU 

-38.I89 

1-8576 

15.676 


Tg - 10 “3 ia* 

-1.5778 

' 6-7238 

-12.093 

0.55690 

-I.8576 

- 0.2277 

13.781 


Ug ■ 1 X )~3 in. 

11.660 

-78.11a 

0.55690 

-8I.510 

15-876 

- 13-781 

19.97* 


'■q * 10 “3 radian 


-38.I89 

-1.^76 

I5-876 

43^-37 

0 

-77.623 

-3B.I89 

T c m 10 “3 in- 


1.8576 

- 2.0077 

-13-781 

0 

- 9-9231 

0 

-1.8576 

Uq ■ 10 "3 iii. 


15.876 

13-781 

* 9 - 97 * 

- 77-603 

0 

-300.31 

*5-876 

ij) * 10 "3 radian 



II 


-38.I89. 

-1.8576 

15.876 

430-37 

Tjj - 10 “3 jn. 





1-8576 

-e.2277 

-13.781 

0 

xij) « 10 w 3 in. 

i 




15-876 

13.761 

* 9-971 

-77.603 

■ 10 ”3 radian 







1 

-38.1^89 

TR = 10 “3 In. 








1.8376 

Uj - 10 "3 Jn. 








*5-876 

Tj - 10 “3 in. 









▼q. *■ 10 ”3 jn. 









■ Ifaroea and 

% 

T u 


% ** 

*r 

% 

Operairfm 

(u») 

(id) 

(to.-l 

1) ( 31 ) (11) 

(11) 

(11) 


t a = 10-3 in. 

Vg - 30 "3 radian 

t b - H)-3 to. 
ng - 10 “3 to- 
Xg - 10 "3 radian 

I.8576 

*5.876 

Tg - 30 ”3 to- 

-2-0277 

13.781 

Ug - 10“3 in. 

-13.781 

*9-97* 

Tp k ID _ 3 radian 

0 

-77-623 

Tp - 30 '3 to. 

-9.9031 

0 

Up - 10”8 to- 

0 

-100. 31 

xjj - 10*3 radian 

-1.8576 

15.876 

Tj » 10 ”3 to. 

-2.2277 

-13.781 

m iO“3 to. 

13-7®- 

19-97* 

Tf » 10“3 to. 
Tg » 10*3 ill. 




- 3 B.M 9 

-1.8576 

*5-876 

- 6 * 9 - 2 * 

-18.056 

-©7.080 

-5.9020 

16.025 


1.8576 

-2.2277 

•0-3-7® 

-1B.056 

-53-95T 

40-333 

-+• 771 ® 

1-3355 


I5.&76 
13.781 
*9-97+ 
-67.080 
->•0. 533 
-126-37 
-1+.66B 
-1-3355 


-5.9000 

-I.5778 

-11.662 

-7-1310 

0 


16.005 

1.3355 

-1-3355 

0 

- 1.8886 





# 


TABES 12 .- GBOUP OEERATIOBB XHf GEROmG-EJSXT METHOD JOE SBGMEHT ABC 


Hg = 454.34»g + 6.723®^ - 78.4jJL.Ug - (R.H.S. in Hg equation) => 0" 

% = 6.7238^ - 12.0937 B + 0.55690^ - (E.H.S. in % equation) = 0 

Tg = -78. 4llyg + O.5569GVB - 84.5ID15 - (E.H.S. In Tg equation) = 0 


Group 

I 

II 

m 

IT 

“''■''^Disp lac ement 

v A n 10 "3 in. 

Vq = 10 "3 radian 

vq = 10 "3 in. 

U£> = 10"3 in. 

Operatlon''''\. 

0 

11 

? 

11 

11 

S’ 

V A “ T c - = 0 

t a = v c = uc * 0 

t a = ¥ C = T C = 0 

(-lo 3 ) x right-hand, 
side in equation 
for: Ifg^ in* -lb 

5.9020 

-38.489 

I.8576 

45.876 

%, lb 

4.5778 

-I.8576 

-2.2277 

13.781 

Tg, lb 

14-662 

45.876 

-13.781 

49.974 

(lO 3 ) x displacements 
of Joint B: 

Vg ; radian 

- 0.026434 

- 0.21631 

0.035554 

O.OI7847 

V B.» 1x1 ' 

-0. 38424 

-O.2397I 

- 0.17353 

1.1763 

v%, in. 

0.19549 

0.74197 

-0.19720 

O.58253 

Besultant 





forces and momenta 
l«b 

-2.6618 

IO.699 

-I.8871 

3.2615 

Hq, In. -lh 

IO.699 

-389.56 

-IO.O93 

- 53-771 

Sq, 

-I.8&71 

-IO.O93 

- 6.7529 

-10.615 

T C , lb 
Up, in. -lb 

3.2615 

0 

- 53.771 

-38.489 

-10.615 

-I.8576 

- 54.199 

45*876 

lb 

0 

I.8576 

-2.2277 

- 13.781 

t d , lb 

0 

45.876 

13.781 

49.974 


o\ 

-3 
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TABLE 13.- GROUP OPERATIONS TS GROWING-ONIT METHOD FOR SEGMENT ABCD 





*c 

BC 

T C 

- -389*56111 - 1Q.Q93^ttt “ 53 * 771 *nr " 
= -10.093a:j-|- - 6<7^9 i ttt “ 10.6l5xjy - 
= - 53 * 771^11 - IO.6151UJ - 5^*1 99*17 ~ 

(R.H. 5 . in Hq) — 0 
(R.H.S. in Rq) = 0 
(B-H.B. in To) “ 0 


Group 

H 

V 

71 

vn 

vm 

'"''^'^Diaplac ement 

(I) - 1 

vj) = 10 "3 radian 

Tj) = 10 '3 in. 

uj) = 10~3 in. 

operation 


» 

11 

t f 
11 

0 

(I) - Tjj = Ujj =» 0 

(I) = wjj * up *= 0 

(I) = Vp . t d = 0 

(-1) x right-hand 
side in equation 
for: He, in. -Xb 

10.699 

-38.489 

I.8576 

45*876 

Be, U> 


-I.887I 

-1.8576 

-2.2277 

13.781 

T C/ lb 


3.2615 

45-876 

-13.781 

49*974 

Magnitudes of (H 
(IH), and (IV 

x n 

I? 

0.021L97 

-0.25291 

■ - 

0.046327 

-O.OO99OI9 

x m 


- 0*53843 

-2.3436 

0.1D521 

O.85346 

X IV 


0.14430 

1.5564 

- 0.32084 

0.76472 

Forces and moments: 
R a , lb 

-0.94505 

6.7928 

-O.74924 

O.77749 

% in«- 

lb 

6.7928 

- 346.88 

- 16.697 

- 43 * 74-5 

% » 


-0.74924 

-15.697 

-5.65OO 

-12.458 

Td, lb 


0*77749 

43 . 7^5 

- 12.458 

-50.817 

in-- 

lb 

0 

-38.489 

-18.576 

45.876 

Rg, lb 


0 

1.8576 

- 2.2277 

-13*781 

Te, lb 


0 

45.876 

13.781 

_ 

49.974 




l 


* 
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IABUE 14. - GROUP OPERATIONS IN GROWING -UNIT METHOD FOR SEGMENT ABODE 


Ng = -346.88xy-j- - 16.69 ^vii ~ ^3*7^5^111 - (E.H*S« In Ng^ — 0 

Kg = -I 6 . 69717 X “ 5- 6500xy-[i - 12 . 458xyj j j - (R.H.S. In I^jj) = 0 

Tg = -4-3 • 7 45^yj “ 12 • ^-j^fayTT - 50 * 8l7^v TTT ~ (R.H.S. In Tg) = 0 


Group 

U 

I 

H 

m 

'---^Displacement 

(V) = 1 

Vg 55 30 "3 radian 

Tg = 10 “3 in. 

ug a 10 "3 in,. 

Op erat ion~''\^_ 

13 T E 53 “ 0 

0 

i! 

ii 

iP 

1! 

S 

00 = vs = % = ° 

(V) = wg = t e = 0 

(- 1 ) X right-hand, 
aide in equation 
for: Ng, in. -lb 

6.7928 

-38.489 

1.8576 

45*876 

%; ^ 

- 0.74924 

-1.8576 

-2.2277 

I3.78I 

Tg, lb 

O.77749 

45.876 

-13-781 

49.974 

Magnitudes of VI, 
VXE, and Till: 

*71 

0. 028042 

-0.16186 

0.037474 

- 0.0022427 

z vn 

- 0.42660 

4.6760 

0.35713 

0.59436 

I YHI 

0.095744 

2.1885 

-0.39100 

O.83963 

forces and moments 





e A» » 

-O.3605O 

4.1055 

-0.31703 

0.19226 

Ng, in. -lb 

4.1055 

-533*92 

-38.099 

- 29.579 

Rg, lb 

-O.3I703 

-38.099 

-49.295 

- 53*^32 

t e , 11 

0.19226 

-29.579 

-53.432 

- 76.322 

%, 11 

0 

-5.9020 

-4.5778 

- 14.662 

11 

0 

16.025 

1-3355 

- 1.3355 



haeis I?.- maaawaBM or total mHFuemanB nt ceohkwet weed 












































ft* 


table 1 6.- caaoup qpmfiobb nr helaiatiqn meedd ids sbmbt abcd 

\ = -7.1310t a + 5.902015 - 4.5778^ + 14.6621^ - (E.H.S. In B -&) = o" 

% c: 5.90£0 t a - 454.34^ + 6-723^ - ' 7 8 *411 x$ - (R.H.S. in %) = 0 

Eg = -4.5779 t a + 6.723815 - 12.093*5 + 0.5569015 - 0R.H.S. in %) = 0 

T b = 14.662 t a - 78.411^ + 0.55690*2 - 84.5101^ - (r.H.S. in Ig) - 0 


Gkroup 

XV 

XVI 

xra 

DiBplaceoBnt 

Vg » 10 ”3 radian 

v c » 10‘3 in. 

Uq *> ID -8 in. 

Operation 

t C" u C’ v D” t D’* u I)” 0 

V C“ U C' ,V D“ T D“ U D“° 

*C “ T C “ YD “ t B = U D = 0 

(-103) x right -hand 
side in equation for: 




r a, b> 

0 

0 

0 

In* -lb 

-38.489 

I.8576 

45.876 

%, it 

-I.8576 

-2.2277 

13.781 

Tb, it 

45.876 

-13-7&1 

49.974 

(lo3) x displacements of 
Joints A and B: 


-0.70895 


y Aj 

4-0195 

1.2253 

Wgj radian 

-0-32255 

0.054293 

-0.014540 

r s , in. 

-1.784£ 

O.O98882 

0.70554 

^ in. 

I.5277 

-0-33579 

0.82205 

Torces and moments: 



-4o.66e 

Hq, in. -lb 

-346.56 

-17.678 

Be, It 

-A7.678 

-5-4150 

-12.928 

T C , It 

-40.662 

-12-928 

-50.203 

Nd, in- -It 

-38.489 

-I.8576 

45.876 

Bp, lb 

1.8576 

-2.2277 

-13-781 

Td, lb 

45.876 

13.7S1 

49.974 
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TABU 17.- (HODP OBERAm® JB BJUMTiail HSEBDD JOB BBSCBS CH5F-0- 


-J 

ro 


Bj 1 -6 49 - 24 *^ - 10.O56tj - 67.060^ - 5.90£0T y + 16.029 t q - (fe.H. 3 . la Kg) » 0 
Bg - -Ifl.056 ^ - ?3-99T»b - JiO-533^ - + i-3355v 0 - (R.H.s. la %) - 0 
■Eg - -^.oeong - H0.533t e - 126.37115 - 14.662^ - 1.3355t 0 - (b.H.B. in Eg) = 0 
% » -5.g020wg - i..577ftr B - I4.66i2i£ - 7.13lOrji - (B.S.8. In %) = 0 
% - 16.02?^. ♦ 1.3353*® - 1- 3335115 - 1.8886r a - (B-H.b. ln%).0 


Grot® 

xvm 

nz 

xz 


Vp m ID “3 radian 

vjj = 10' 3 la- 

lip « 10 “ 3 la. 

Operatic* ! 

v c = T C = ^ “ T D - ” 0 

y C“ r C mV C my T> mU D mi:> 

lr C- T C“ 1J c =-si ) =ii D = 0 

(-lt> 3 ) x rl^it-hand 

side in oqpatlon: 

-38.489 



% in. -lb 

-1.8576 

45.676 

Hg, It 

1.8576 

-2.2277 

-13.781 

Tfc It 

45.876 

13.781 

49.974 

%, It 

0 

0 

0 

'* 0 , It 

0 

0 

0 

( 1 £> 3 ) x dlsplatcemanta 
of Joints I, /, and G-‘ 



-O.OO8775 

*g, radian 

-0.1703 

-0.032775 

T J» la- 

-0.4297U 

-0.19092 

-O.78327 

Hj, la. 

0.74456 

0.23168 

0.79421 

Tj, la. 

-1.11^1 

-0.32667 

-1-1229 

▼a» 

- 2.2753 

-0.57676 

-1.1900 

Jorcea sod noaects: 




ir c , la. -lb 

-38.489 

1.6776 

45.676 

xt> 

-1.8776 

-2,2277 

13.781 

V 111 

J15.876 

- 13.871 

49.974 

*n, la. -lb 

-392.46 

u.535 

- 42.305 

it 

11.535 

- 6.2441 

12.706 


-42.305 

12.706 

-50.259 
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TABLE l8. - GROUP -OPERATIONS TABLE FDR RELAXATION METHOD 
[Forces and. moments at Joints A f B, E # F, and G are zero for all operations^ 


^sForces and 
x s mcments 

OperatlcnN. 

% 

(in. -lb) 

*C 

(it) 

T c 

(ID) 

*D 

(in. -lb) 

e d 

(ID) 

(ID) 

(XV) = 1 

-346.45 

-17.678 

-40.662 

-38.489 

I.8576 

45.876 

(XVI) = 1 

-17.678 

-5.4150 

-12.928 

-I.8576 

-2.2277 

13.781 

(xvii) * 1 

-40.662 

-12.928 

-50.203 

45.&76 

-13.781 

49.974 

(xvm) = i 

-38.489 

-1.8576 

4 5 .Q 76 

-392.46 

11*535 

-42.305 

(XXX) = 1 

1.8576 

-2.2277 

-13.781 

H-535 

-6.2441 

12.706 

(XX) = 1 

45.876 

13.781 

49.974 

-42*305 

12.706 

-50.259 













TABU 19*- ISLAZXEnV 



Choafc tabl« 0.001 0-310 O-OSg 


St»p 

acd 

Operation 

■b 

pi) 

T » 

(lb) 


1 

tetomoi icodi 
302 x (I) 

- 9 « 

ms 

-2» 

326 


2 

-lfl.lt. X (M) 

052 

a 

0T 

0 


3 

Zllt X (K) 

Ql 

67 

-El 


V 

1(.08 X (0) 

8 

55 

1 (6 
0 


? 

■J- 3 -T X (K) 

63 

-63 

« 

-17 


6 

-0.76 x (xm) 

0 

uo 

5 


ebook tablo 

10 

3.390 

-9 , 

-10.650 
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Figure 1.- Reinforced panel with conditions at both ends specified in 

terms of force. 



Figure 2.- Reinforced panel with conditions at one end specified in terms 
of force and at the other in terms of displacements. 
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Figure 6.- Circular ring with antisymmetric loads. 

500 lb 



Figure 7 .- Egg-shaped ring with symmetric loads. 
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Figure 8 . - Oval -shaped ring with positive directions of forces and 

moments shown. 
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